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Longitudinally polarized, unpolarized and forward-backward mass differential cross sec- 
tions for Drell-Yan lepton-pair production by arbitrary vector bosons are calculated in next- 
to-leading order (NLO) QCD. Analytical results are presented in a form valid for all consistent 
n-dimensional regularization schemes, with the mass factorization scheme kept general. NLO 
predictions for all Drell-Yan type processes (W , Z and 7*) at BNL's relativistic heavy ion 
collider (RHIC) are made using polarized parton distributions which fit the recent deep-inelastic 
scattering data. These are examined as tools in the determination of the polarized parton dis- 
tributions and the unpolarized u/d ratio. NLO predictions for the forward-backward lepton 
asymmetry at Fermilab are made and the precision determination of sin 2 8w from future runs 
is studied. In all the above, the QCD corrections are found to be significant. An introductory 
discussion is given of various theoretical issues, such as allowable factorization and regularization 
schemes, and scale dependences. 
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I. INTRODUCTION 

1 Two major areas of interest within the standard model are the determination of the polarized parton distribu- 

tions of the proton and higher precision determinations of the electroweak mixing angle, sin 2 9w , as a constraint 
on the Higgs boson mass and new physics. One process useful in exploring both areas is Drell-Yan lepton-pair 
production. We shall present a clear picture of the Drell-Yan process at one-loop in QCD within a general frame- 
work which should be both instructive and useful. We will do so by considering the general interference between 
two vector bosons with arbitrary mass, width and couplings, which decay into a general lepton- antilepton pair 
(including neutrinos) and are produced via quark-antiquark fusion (one of the (anti)quarks may of course arise 
■ from an initial state gluon). In this way, we may consider possible new physics contributions, such as Z' bosons 
and four-fermion interactions, by appropriate choice of couplings, etc. . . 

The emphasis here will be on presenting complete analytical results in a form valid for all consistent n- 
dimensional regularization schemes and within a general mass factorization framework. In addition, we will 
Qh. consider all possible longitudinal polarization states of the initial hadrons. Mass differential cross sections and 
asymmetries will be presented and the effect of the next-to- leading order (NLO) subprocesses will be highlighted 
and explained pedagogically. We will also discuss various constraints on allowed regularization and factorization 
schemes and describe the origin of the scale dependence of the one-loop corrected predictions in a general fashion. 

The one-loop QCD corrections to the longitudinally polarized Drell-Yan process have been studied in several 
papers Q|- ^| using various regularization prescriptions (for the 75) and factorization schemes while considering 
the production by specific bosons. Here, we keep the formalism completely general. Results are kept in a 
simple form by considering explicitly only mass differential cross sections. As an application, we will study how 
lepton-pair production by 7*, Z and bosons can be used to extract the polarized parton distributions from 
forthcoming planned polarized pp collision experiments at BNL's relativistic heavy ion collider (RHIC), which is 
scheduled to start running in 1999. Longitudinally polarized pp collisions with both beams polarized are expected 
to begin in 2000 Q with sufficient muon coverage to perform precision spin studies, in particular using W^'s. 

The QCD corrections are necessary in order to reduce the process dependence of the parton distributions. 
This allows comparison with and use of those distributions obtained in polarized deep-inelastic scattering (DIS). 
This form of global analysis will prove invaluable since, with the Drell-Yan process, we are sensitive to the 
sea-quark distributions, which are currently almost totally unconstrained from DIS. We can make use of the 
DIS determinations of the polarized valence distributions however, since we work consistently at NLO in QCD. 
Eventually RHIC will be able to improve those valence determinations due to the increased flavor separation in 
production. 

Another use of Drell-Yan at RHIC is a precision determination of the unpolarized u/d ratio at fairly large x and 
at very high energy scales where other experiments fall short. This can provide information on the limit x — » 1 
at lower energy scales, via perturbative evolution. Similar statements apply to the polarized parton distributions 
with regard to the large- X sensitivity. All predictions will lie within the energy range 100 GeV < ^/s pp < 700 GeV, 
although RHIC, as currently envisioned, may not be able to run much above 500 GeV. 

The Drell-Yan process also provides a useful way of obtaining s'm 2 9\y- One may ask why we would try to 
measure sin 2 6w in hadron-hadron collisions when good measurements exist in Z production from e + e~ anni- 
hilation at LEP and SLAG The main reason is that these measurements are nearly complete. Future SLAC 
measurements will improve the earlier SLAC value, but a large discrepancy with the LEP measurement presently 
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exists and should not be expected to vanish. Thus, what is needed is an independent high precision measure of 
sin 2 9w- With various high luminosity scenarios planned for Fermilab's Run II |7|, such a high precision mea- 
surement is indeed possible. This is because the forward-backward lepton asymmetry (in the lepton-pair rest 
frame) is sensitive to sin 2 8w In a best case scenario, one could surpass the present (and future) average from 
LEP and SLAC combined. In a worst case scenario, a measurement at the same level of precision as the present 
SLAC measurement should be attainable. Various studies HQ| have been done to show how such a high precision 
determination would significantly constrain the allowed mass region for the standard model Higgs and provide a 
probe of new physics. 

QCD corrections to the charge asymmetry of lepton pairs produced in pp collisions have previously been 
investigated H . Unfortunately, the physical observables considered are not the ones used in determining sin 2 9w 
and there does not appear to be any way to straightforwardly convert them into a useful form. Very recently, 
a paper appeared performing detailed numerical studies of the QED and QCD corrections to the forward- 
backward Drell-Yan asymmetry in pp and pp collisions. There, Monte Carlo methods were used and similar 
results were obtained, even though a different definition of the forward-backward asymmetry beyond leading 
order was used, apparently in order to minimize the QCD corrections. As well, recently, soft gluon resummation 
effects on the lepton angular distribution from the decay of Z's produced at the Fermilab Tevatron were briefly 
considered in [jLl| . Here, complete analytical results are presented for the mass differential cross sections and 
these do not appear elsewhere to our best knowledge. Also, hadron polarization effects as well as regularization 
and factorization scheme dependences are presented explicitly, which is not done elsewhere. One also sees quite 
clearly the structure of the QCD corrections and the origin of that structure. 

The paper is organized as follows. In Section II we present our general formalism and describe the observ- 
ables being considered. In Section III we discuss various features of dimensional regularization and dimensional 
reduction. In Sections IV-VI we compute the (singular) subprocess cross sections. In Section VII we perform 
the factorization of the mass singularities and discuss constraints on allowable regularization and factorization 
schemes. We also discuss how to convert subprocess cross sections and parton distributions from one scheme 
to another. In Section VIII we present the final analytical results and discuss the scale dependence issue in a 
general fashion. In Sections IX-XI we present numerical results relevant to RHIC and examine the sensitivity to 
the polarized parton distributions as well as the unpolarized u/d ratio. In Section XII we present the forward- 
backward Drell-Yan asymmetry relevant for the Fermilab Tevatron and discuss the extraction of sin 2 9w from its 
measurement. Finally, in Section XIII we present our conclusions and summarize the work. Throughout, we have 
tried to avoid giving standard discussions of general issues in order to present new and more general perspectives 
which should be useful to non-experts. 

II. GENERAL PROCESS AND FORMALISM 

The Drell-Yan process with initial hadrons A, B of definite chirality is 

A(P U X A ) + B(P 2 ,X B ) -» Kp 3 ) + T(pa) + X, (1) 

where Xa, Xb denote chiralities, 11 represents a lepton pair and X is an arbitrary hadronic final state. We have 
the following invariant observables, 

S=(P 1 +P 2 ) 2 , M 2 = q 2 , t = M 2 /S, (2) 

where 

Q = P3+P4- (3) 

Denoting the y-axis as the direction of motion of hadron A, we may also define, in the cm. of A and B, the 
observable 
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The general 2^2 [2 — > 3] parton subprocess contributing to (Q) may be written as 

a(pi, Ai) + b(p 2 , A 2 ) 

- B*(q) + [c(k)] -» l(p 3 ) + l( Pi ) + [c(k)], (5) 

where B{ = 7, Z, and I = l~, v\\ I = l + , vi (in the standard model). We must consider the cases a = q, b = q, 

(-) (-) (-) (-) 

c = g; a = q, = q, c = g; a = q , b = g , c = q ; a = g , = q , c = q . 

Using the parton model relations 

pi = x a Pi, p 2 = x b P 2 , (6) 
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we may define the subprocess invariants, 

x2 a M 2 M 2 T 

s=(p 1 +p 2 ) = x a x b S, w= = - = , (7) 

S oX a X b X a X b 

where we took all external momenta to be massless, as usual. The parton momentum distributions are given by 

Fjf I (x i ,fi 2 )=Xifl /I (x i ,fj, 2 ), k = u,l, (8) 

where 

fj'ixun 2 ) = ff /+ {x u M 2 ) + fH; + ( Xi ,^) (9) 

and 

fi'^xuf) = /; 7 ; + (^m 2 ) - rlj+^n 2 )- (io) 

Here f^/x^x-i, p 2 ) is the probability of finding parton i with chirality Xi and momentum fraction Xi in hadron I 
having chirality Aj, evaluated at renormalization scale p 2 . 
Introduce 

z = c 0S r, (ii) 

where 9* is the angle between p 3 and Pi (or, equivalently, pi) in the 11 rest frame. Then we may define 
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1 da f , da ,„ . 

dz T z ± J_ 1 dz T z - (12) 



Hence a F+B has the interpretation as the usual leptonic integrated cross section. We call a F ~ B the forward- 
backward cross section. 

Let a ab (X a , X b ) denote the cross section for colliding partons (or hadrons with some degree of polarization) a, b 
having chiralities (or beam polarization directions) A , A;,. Then there are only 4 combinations of A a , At, leading 
to well factorized parton model expressions, as may be straightforwardly verified. All other observables may be 
expressed in terms of these. They are 
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The notation is straightforward: a^ n denotes the cross section when a has polarization m and b has polarization 



u = (A = +) + (A = -), i = (A = +)-(A = -), (17) 

so that u denotes unpolarized and / denotes longitudinally polarized. The factor 1/4 is required so that <r uu has 
the interpretation as the spin averaged cross section. 

The parton model expression for the mass differential Drell-Yan cross section, for general beam polarization, is 

r ]„AB,F±B pi pi j-ab,F±B 

MO .mm ~ A ~ R X > I / , „r, I A / 0\ «/. / R / 0\ MO ™„ 



(IM S^S^JdXaJ d Xb fZl A {x a ,H 2 )f b J B {x b ,H 2 )- 



SiS^f—t -K /A (^^)F b / B (x b ,p 2 )^^ (18) 
SiS B Y, f T dXF0 f ^--^—F^(x a ^ 2 )F b J B (x bl ^- 

ab J-(1-t) Jw! w x a T X b 



ab,F±B 



where 



r x /x 2 fq +4t/w±x F0 

X F0 =X a -X b , W! = - : r, X a . b = . (19) 

1 - \xfo\ * 



a ab is the subprocess cross section corresponding to (||) and 

Si = 1, Si = V 1 , (20) 

where V 1 is the degree of polarization of the beam of hadrons /. The only assumption made in the derivation 
of ( |l8| ) is that the degree of beam polarization is the same in both polarization directions. The same factors 
of beam polarization will always enter, regardless of the specific differential cross section being considered. The 
experimental results (for the a mn ) will have to be divided by the appropriate factors of beam polarization in 
constructing asymmetries. In general, this is unavoidable since the beam polarization will vary as a function of 
time for any given experiment. In numerical calculations taking evenly spaced integration points (i.e. Simpson's 



rule) the third form of (18) is most convenient because x a and Xb are integrated over symmetrically while the 
inner integral is over w and da/dM contains functions singular in 1 — w (i.e. 5(1 — w), 1/(1 — u>)+, . . . ). 

We will consider lepton-pair production by N arbitrary vector bosons, having masses Mi, widths Ti and 
couplings to fermions given by 

"•;-"(</:,/• fl :,,-;). (21) 

Hence we may write 

dM dM ' [ > 

where ff'^g is the total interference contribution from Bi, Bj (i.e. for i ^ j it is the sum over the i,j and j, i 
contributions). For M, T, — > we obtain the photon contribution so that the argumentation is completely 
general. This formalism applies to Z and W production and the Drell-Yan production of any bosons whose 
coupling to fermions is given by ( pl| ) , provided there are no new contributing diagrams. This includes non standard 
model Z' and W^' contributions. As well, we can include four-fermion interactions. For production, the A 
will depend on both quark flavors at the vertex via the appropriate CKM matrix element. 

The normalization convention of our squared amplitudes is easily inferred from the standard leading order 
parton model cross section expression, 
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dMdxpdz irM 3 x a + Xb 



^£^/ A (w 2 )^ /s ( 



X b ,fi 



><s2\M\ 2 m t BiB ., (23) 

i>j 



where 



ab = qq, qq, x F = x F0 (24) 

and BiB . denotes the net interference between the amplitude involving Bi and that involving Bj. This 

corresponds to the convention ^ A u(p, X)u(p, A) =^/2 for p 2 = 0, which we maintain for consistency with Q. 

As far as QCD corrections are concerned, we will only consider da/dM. In imposing exerimental cuts, one 
restricts the phase space of the outgoing leptons. Since the QCD corrections effect the hadronic sector, the 
features of the corrections should not be greatly influenced by the leptonic cuts. In the limit of small cuts (i.e. 
good leptonic coverage) the expressions in this paper become exact. We will not address further the issue of 
leptonic cuts in what follows. 



III. n-DIMENSIONAL REGULARIZATION SCHEME PECULIARITIES 



Since we will need to present the leading order (LO) cross sections in n (= 4 — 2e) dimensions in the next 
section, it is appropriate to discuss the various n-dimensional regularization schemes at this point. We will 
consider regularization by dimensional regularization (DREG) |12]]l3]| and by dimensional reduction (DRED) 

In |]J , the details concerning DREG and DRED were summarized and will not be repeated here. The important 
points will be discussed however. Within DREG, there are two commonly used schemes for dealing with the 75 
matrix (and £^9° tensor) which arises in polarized processes. In the 't Hooft-Veltman-Breitenlohner-Maison 
(HVBM) scheme |l^Jl^], all quantities are mathematically well-defined, but the non-anticommuting 75 leads 
to physical problems. These necessitate finite renormalizations of the polarized parton distributions and UV 
counterterms (for non QCD/QED vertices) as will be discussed throughout this paper. In the anticommuting-75 
scheme [ fL6| , massless quark helicity is conserved, but only at the expense of mathematical consistency when an 
odd number of 75 's arise in the traces. Also, one needs to devise prescriptions jl7j for dealing with the e^ v P a 
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tensor (arising from polarized gluons), since nothing other than the HVBM definition has been proven to be 
mathematically consistent in DREG. 

Even when performing calculations where only an even number of 75 's occur in the traces, one may wonder if 
the result is meaningful in general, using an anticommuting 75 . For the case of longitudinal polarization, the 75 is 
connected with helicity, and we understand the need for helicity conservation on physical grounds. For processes 
such as the transverse Drell-Yan process, the 75 also arises. Here though, there is no physical motivation for 
using an anticommuting 75 in DREG, other than that it simplifies calculations. More study is required in this 
area since recently two-loop transversity splitting functions have been calculated by three groups ]l8|-p0| using 
an anticommuting 75. The transverse case is somewhat more subtle than the longitudinal case, however, due to 
the additional axis which enters. We do not consider transversity further in this work. 

So, in general, the anticommuting-75 scheme |]l6| should only be used in situations where it is physically 
motivated. And in those situations, it should always be shown that the same results can be obtained using the 
HVBM scheme via finite renormalizations of the parton distributions or by the addition of UV counterterms. The 
above renormalizations and UV counterterms should be unique (in any given gauge) for some class of subprocesses. 
This is the approach we will follow in this paper. All results are understood to refer to consistent schemes: those 
schemes which can be related to the HVBM scheme as described above. Hence, we will always refer to the HVBM 
result as being the DREG result. 

More recently, a 75 scheme was introduced |2lJ which reproduces many of the desirable features of the HVBM 
scheme, such as the correct ABJ anomaly, and claims mathematical consistency without violating helicity con- 
servation. This scheme uses a reading point (a specific leftmost 7-matrix in the trace) and maintains an anticom- 
muting 75, which requires non-cyclicity of the traces. What remains, therefore, is to show that this scheme leads 
to process independent n-dimensional splitting functions. This will be discussed below. 

Because of the tedious nature of the HVBM scheme, and the other problems mentioned earlier, it is more 
straightforward to first use DRED to calculate subprocess cross sections, then present them in a form valid for 
all consistent n-dimensional regularization schemes using the technique of [Q . The only drawback is that one has 
to add a UV counterterm to the quark-7(Z, W^) vertex when using DRED. Fortunately, this counterterm is well 
established and unambiguous. 

The DRED result for a physical cross section is (formally speaking) defined as the result obtained by contracting 
all tensors and taking all traces in 4 dimensions, then performing all phase space integrals in n dimensions (and 
adding any necessary UV counterterms). When considering QCD/QED corrections, one needn't actually take 
the traces first. One may work directly at the amplitude level for loop diagrams. For other theories, however, like 
supersymmetry, the ill-definedness of g £ v for noninteger n < 4 leads to an ambiguity such that the results may 
depend on the order of operations 0] . In the HVBM scheme (n > 4) , is defined since there are an infinite 
number of integer dimensions with n > 4. Hence we may continue g^ v to n > 4 by simply using relations valid 
in integer dimensions, but take the continuous limit n — > 4 on the real axis (or complex plane) at the end of the 
calculation (i.e. after all divergences have been cancelled), as if one compressed the infinite number of dimensions 
into a continuum for n > 4. In practice, of course, one continues to noninteger n from the beginning rather 
than at the end. As there are a finite number of dimensions with n < 4, we cannot continue g^ v into the region 
n < 4. Hence, strictly speaking, we must first perform all contractions and traces in 4 dimensions, then go to 
n dimensions to perform the phase space integrals, since only scalar products of vectors will remain, which are 
manifestly continuable to n dimensions. In that way, g£" will never be contracted into traces containing 75 or 
with s lLV( " T , which is where the ambiguity could arise. 

In order to present da^^ ±B jdM in a form valid for all consistent n-dimensional schemes, we must first give 
the general form of the n-dimensional Altarelli-Parisi |g3| splitting functions, PJj(z), related to the probability of 
parton j splitting into a collinear parton i having momentum fraction z, plus an arbitrary final state carrying the 
rest of the momentum. The reason is that the regularization scheme dependences arising from mass singularities 
are entirely contained in these functions, as was shown in Q (and references therein). Hence, we discuss them 
here. We may write 

P%(z,e) = lf j (z)+ePe j (z), (25) 
where Pij(z) is the usual 4-dimensional splitting function. In DRED, 

P£ DRED (z,e) = i*(z). (26) 

We will also use the notation, 

P$(z, e) = P%<(z, s) + 5(1 z)P% s (e), (27) 

and define 



A u Py = Py j ++ + Pij, — 1_, A/Py = Py j++ — Pij ! — 1_, (28) 



where the +, — denote the respective chiralities. 

In DREG, the unpolarized one-loop splitting functions are given by [E4J 
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A u P*(z) = C F 



(l-z)+ 



- 1 - z 



A u PUz)=C F 



-(!-*) + -*(!-*) 



i(l-2z + 2z 2 ), A„P| fl (z) 
C F 



(29) 



H-2. 

z 



^uP £ qq (z) = -C F z, 



= 2N, 



c 



1 



(!-*)_ 



1 

- -2 

z 



z(l 



^6(1- z), A u P* g (z) = ^6(l-z), 



where (3o = -y^c — %N F and the usual convention of 2 — 2e (= n — 2) gluon polarization states was used. In 
fact, the latter is somewhat more than a convention, since it is justified on physical grounds. 

Within DREG, we need only consider longitudinally polarized splitting functions determined in the HVBM 
scheme since, at present, it is the only established mathematically consistent scheme therein, which also leads to 
process independent n-dimensional splitting functions. At one- loop they are given by 



A,i* (*) = A u P± Q (z), A,i* (z) = C 



89 



3(i-*) + ia(i-*) 



^P qg (z) 



1/2, A,i* (z) = -(1 - z), 



(30) 



A,P„ 4 „(z) = C F (2 - z), A^Jz) = 2C F (1 - z), 



A,i* (z) = 27V C 



(l-z)_ 



- 2z+ 1 



^(1 - z), A ; P; g (z) = 47V C (1 - z) + ^5(1 - z), 



as can be inferred from |25| . 

These n-dimensional splitting functions were determined by investigating the factorization properties of the 
NLO squared amplitudes in various collinear limits for various subprocesses. In pq] it was found that one always 
obtains a Born term multiplied by the appropriate n-dimensional splitting function in the limit where two partons 
are collinear, leading to a pole in a propagator. It is important that the Born term be exactly the n-dimensional 
one appropriate to the regularization scheme in question; in the above case, HVBM. We investigated what happens 
in the scheme of for the qg subprocess of longitudinally polarized Drell-Yan. For simplicity, we chose the 75 
coming from the polarized incoming quark as reading point. We found that the wrong Born term arose, in the 
collinear limit where the outgoing quark is collinear with the initial state gluon. Namely, that of HVBM rather 
than that of p]f| . The basic reason is that the squared amplitude for the qg subprocess in the scheme of j2]J is 
equivalent to the HVBM result, with our choice of reading point. Hence, in the collinear limit, the Born term 
(qq) and n-dimensional splitting function (AiP^ g ) which arise correspond to the HVBM ones. Unfortunately, 
the Born term for the qq subprocess in the scheme of pl[ demands an anticommuting 75, hence the result is 
simply minus the unpolarizcd one. This is not true in HVBM, due to the non-anticommuting 75. Hence the Born 
terms are different in the two schemes and the behaviour of the qg squared amplitude in the collinear limit is 
unphysical. This physical inconsistency is clearly unacceptable and it will lead to process dependent n-dimensional 
splitting functions. A different choice of reading point might rectify the situation, but such an ambiguity is also 
unacceptable. Unfortunately, the prescription of pl| gives no unambiguous definition of what the reading point 
should be in all cases. Therefore some extension is in order before it can be applied in practical QCD corrections 
to polarized processes. We will not discuss the scheme of ]2l| in further detail. 

The only subtlety associated with calculating the AjP,™ in the HVBM scheme is that one must additionally 

perform an integration over the k momenta (i.e. the components between 4 and n dimensions), rather than simply 
taking a collinear limit [Q] as is done in anticommuting-75 schemes and in the unpolarized case. One still obtains 
process independent splitting functions since the collinear phase space structure is process independent. Even in 
some anticommuting-75 schemes (l7| , one must resort to Sudakov kinematics in the collinear limit, making the 
process nontrivial. In that case, one must drop certain potentially finite terms of 0{e) which only vanish after 
integration. 

In what follows, we shall present tensors and squared amplitudes calculated using DRED. Then, using the 
method of we present the necessary one- loop subprocess cross sections in a form valid for all consistent n- 
dimensional regularization schemes. By "n-dimensional regularization schemes" , we mean schemes in which all 
divergences are regularized by dimensional continuation. This excludes all forms of off-shell and cutoff regular- 
ization schemes which also involve dimensional continuation to regularize some of the divergences. The mass 
factorization scheme will be kept general. 



IV. LEADING ORDER CROSS SECTIONS 



Unless otherwise stated, the expressions presented here are calculated using DRED. The unintegrated leptonic 
tensor Lg B , is defined as the product of L B and L B , , where L B is the leading order amplitude consisting of a 
lepton pair attached to boson B at a vertex having index a. It is given by 



G 



fl 2e 

T OLD h 1 /[/ I | / \rrsOL3 

l bb> = — c i c i[{gvigvi + gaig a i) T u 

-(gaWvl + 9vl9al)T%f], (31) 



with 



T 2 f = 4*£ Q ^> 3 ^4„ (32) 

where the arbitrary mass scale [i 2e arises from the n-dimensional coupling e 2 — > e 2 /i 2£ . Since [i is arbitrary, the 
physical predictions should not depend on it, order by order in a s . We will show how this is explicitly satisfied 
at one- loop in section VIII. 

Now define the integrated leptonic tensor as 

= f Vp 3 

Jf+b 



ra p _ f d n 1 p 3 5[{q-p 3 ) 2 ] p 

Lf + b ~ J (27T)"-i2p 3 ,o q 2 ' [ ' 

Here we omitted the BB' indices, and will often do so for compactness. One finds |Q 

(34) 



r a/3 



a <x0 <x0 

(l-s) £ ^- + ^--(3-2 £ )^- 



where k now generalizes to 



, M 2£ (g 2 )- £ r(l-g) (g v ig' vl +g a ig' al ) 
C;C 44- 2e n 2-e r(l-2e)(3-2e)(l-2e)' 1 j 



It is interesting to note that the corresponding DREG tensor, in an anticommuting-75 scheme, is obtained by 
replacing g af3 — > g" 13 . This gives 



i -F+S,DREG — "A 1 — £ . 



SLi__ fl c/3 
2 y " 



(36) 



Of course, is not defined in such a scheme, but it does not contribute here. In all schemes, the part ~ q a q 13 
does not contribute to the cross section, as follows from gauge invariance. is thus effectively a constant 

tensor. We nonetheless keep all the terms for completeness. 

We may define the forward-backward integrated leptonic tensor with respect to an arbitrary massless vector, 
p, in a covariant fashion using 



JF+B 



I F+B 

-e(p-P3~p-Q/m aP (37) 

V P3 L a() (38) 

F-B(p) 

Kie a ^ u ^^, (39) 
2p-q 



where 9 is the step function and 



u 2e (a 2 )~ £ Y(l - e) 

« = w'l 05-2, J-, r(1 _ 2e f (e)(9aig' vl + gvig' al ), (40) 



F( £ ) = 4-L— -i/ dyiy-^l-y) 1 ^ 



with 

= A 

r 2 (l-e) 

-y 1 -°(l-y)- B ] = l + 0(e). (41) 

The most general tensor structure of ^p^_B(p) ma y be obtained by noting that only the antisymmetric part of 
L a P contributes. Hence ^-p^B(p) mus ^ be antisymmetric in a, (3 and can only depend on the momenta p and q. 
Then one uses the usual projection methods, as was done to obtain ([mJ). We will use the notation 



JF-B JF-B( Pl ) 

It is easy to check that ([57]) leads to the usual definition of the forward-backward asymmetry, Eq. (|l2|), by 
working in the rest frame of q. Choices of reference axes other than p\ are of course possible, but using p\ allows 
a straightforward, covariant treatment. Having performed the forward-backward leptonic integration analytically 
greatly simplifies the next-to-leading order calculation. 

We may next define the subprocess hadronic tensor W"f BB , through the subprocess squared Feynman amplitude 
contribution from B,B' interference, 

\M\ 2 ab>BB , ee D BB ,L B a fw^ BB ,{2 - S BB ,), (43) 

where 

n -i = (Af 2 - M 2 )(M 2 - M 2 ,) + M B M B i T B T B i 

BB ' ~ [(M 2 - M|)(M2 - M%) + M B M B >T B T B '} 2 + [M B T B (M 2 - M%) - M B ,T B ,(M 2 - Mf )] 2 ' 1 > 

For a = q,b = q, the leading order hadronic tensor is 

W q f = ^c q c' q {{ 9vq g' vq + g aq g' aq )[T^ (1 - A x A a ) + T 2 f(X 1 - A 2 )] 

-(g aq g' vq + g^g' aq )[T^(Xi - a 2 ) + r 2 f (1 - AiA 2 )]}, (45) 

where 

M 2 

T x f = 4(p?pg +p£p? - — 3 Q/3 ), T 2 f = -4ie^»p llM p 2 „. (46) 



Using (|31j) and ( [45|) in ( |43| ) gives the LO result for a = q, b = q, 

\M\ 2 qq . BB ,(z) 



= (2- S BB/ )D- B , ^ c;C;C g Cg{(l - AiA 2 )[(.g l ,zg^ + gaig' a i){gv Q g' vq +9a q g' aq ){l + z 2 ) 

+ 2 (9aig' v l + 9vig'al){9a q g' vq + 9v q g' aq ) z \ + (A 2 - Al)[(^^^ + S'aJfl'li) 

x (9aq9' vq + 9v q g' aq ){l + z 2 ) + 2(g a ig' vl +g v ig' al )(g vq g' vq +g aq g' aq )z}}. (47) 

This reproduces the result of J2(| which includes the effect of Z' bosons in the Drell-Yan process. We may express 
the above in a covariant fashion using 

z = (u 2 - t 2 )/M 2 , l + z 2 = 2(t 2 + u 2 )/M 2 , 

t ee (p! - P3 ) 2 = (p2 -p 4 ) 2 , u= (p 2 -Pa) 2 = (pi -P4) 2 - (48) 



M F±B 

JVl ab.BB> 



Defining 

/ Vp,L B J'w: b % B ,=C B f^ B W% BB ,, (49) 
we have, for a — q,b = q and in leading order, 

M qq +B = 2 3jv c c i c 'g K ( 2 ~ £ )\ig^9'vq + gaqg' aq )( l - ^2) + (g aq g' vq + gv q g' aq )(M - Ai)] (50) 



fJ- 2E M 2 111 11 

M qf B = OS at c 1 c 'qK[{gaqg' + ~ ^2) + (fftiaffL + SagflL )(A 2 " Ai)], (51) 



and 



where we used ( |34| ) and (|39|), respectively, in (|49j). The correctness of this approach (i.e. integrating over the 
leptonic tensor first, then contracting the leptonic and hadronic tensors) was verified by contracting first, then 
performing the leptonic integration. One obtains exactly the same result, in n dimensions. 
This is straightforwardly checked by noting that, in the rest frame of q, 

where 
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d U! 3 



dBi sin 1 "" 26 Ox I d9 2 sm-^6 2 
ii Jo 



(53) 



and the 9i represent the first two of the n — 2 angles of p 3 in n dimensions |27|] and &\ is taken to be the angle 
between p 3 and p\ (i.e. 9*). Since |M| 2 does not depend on 9 2 , 



d lo 3 — > 2tt 



r(l - 2e) 

r 2 (i-e) 7o 

We simply make the substitutions, 

z = y-(l- y), 

in (|47|), then integrate over y. 
The 2 — > 2 phase space is 



dyy E (l-y)~ 



1 + 



1 



%(l-y), 



da 



ab,F±B ' 
,.BB' 



dM 



LO 



X has the form 



n-1 A/f afc - F ± B /o X \ — ab,F±B , U /, 



X m ; Ri3'( £ ) = ( 2 ~ Obb'J^-t? n ^ss' i? ~ 



2tt7V c £>BB' 



Here the K l BB , ,F±B are leptonic factors given by 



K 



rr,F+i 

BB' 



= o(9vig' v l +9aig'al) 



3 (j, Ae tt £ r(l-e) 



rnn.BB' ' 



2-e 



2 2- 2e q 2s r(l - 2e) (3 - 2e)(l - 2e) 



(54) 



(55) 



(56) 
(57) 

(58) 
(59) 



A' 



i + r,F-B 

BB' 



= ^(9al9vl + 9vl9al) 



fi 4s if r(i-g) 

2- 2e ^7r(l - 2s) 



F(e) 



(60) 



The factors in the square brackets are equal to l + (D(s) and hence may be taken as 1 for calculational purposes, as 
they will always factorize exactly in next-to-leading order. Similarly, the Born term may be computed analogously 
in any valid n-dimcnsional regularization scheme. However, the Born term always factors out of the singular parts, 
which cancel (after mass factorization). Hence, the scheme dependence of the Born term cancels in the limit e — > 0. 
This will be demonstrated explicitly in the following sections. 
The F ab are hadronic factors whose expressions are 



qq,F+B _ _^qq,F+B _ pqq,F+B _ _ p qq,F+B _ 



T^qq.r -rr> _ ZT" 

r uu.BB' ~ II. BB 



,BB 



E> IT" 



II. BB' 



9vq9vq 1 9aq9aq' 



(61) 



and 



w qq,F+B _ _ p qq,F+B _ p qq,F+B _ _ p qq,F+B / / 

ul.BB' — r lu.BB' " r lu,BB' ~ ul,BB' — 9aq9vq T 9vq9 aq , 



qq,F-B _ _pqq,F-B _ 



uu.BB' 



F 



qq,F-B 
ul.BB' 



ILBB' 



-f: 



qq,F-B 



lu,BB' 



rpqq.F—B r?qq,F '— B f . i 

~ ~- t UU,BB' ~ ^ll.BB' ~ 9aq9vq ' 9vq9 a qi 



p qq,F~B _ p qq.F-B 
' r lu.BB' — r ul.BB' 



9vq9vq 9aq9aq- 



(62) 

(63) 
(64) 



V. VIRTUAL CORRECTIONS 



The string of gamma matrices, 



(65) 



arises in the calculation of the vertex graph (the massless self-energies vanish in n dimensions). In order to get 
the correct form for the virtual corrections, we must satisfy the relation, 



7 Al 7 p 7 Q (ff t . ? - gaqish^ju = i tl Yi a i' I iu{gvq - gaqis), 

9 



(66) 



which holds for an anticommuting 75. The relation (|66|) is necessitated by the fact that the QED Ward identity 
]2q] between vertex and self-energy graphs holds for the vector part. Therefore, we must be able to anticommute 
the 75 out of the vertex correction so that the axial-vector part factors correctly and the identity is trivially 
satisfied. In non-anticommuting-75 schemes, this relation would be violated by terms of 0(e) (giving rise to a 
finite piece when multiplied by the 1/e UV singularity) which would have to be removed by UV a countertcrm 
analogous to that of DRED which corrects the j-q vertex. Then, the result would be equivalent to that obtained 
using an anticommuting 75 in n dimensions, except that, in general, there could be additional terms of soft origin 
which must cancel when the bremsstrahlung contributions are added. If they did not cancel, the scheme would not 
be physically consistent since soft divergences cannot give rise to scheme dependences in physical cross sections. 

It should be noted that we are assuming a vertex having the form (|2l|). It was shown in |^9| that using a 
symmetrized vertex can avoid the spurious UV terms, in the HVBM scheme, for the vertex b — > s + H. We 
investigated whether such a procedure would work here. Unfortunately, the situation is somewhat different than 
in b — > s + H . Considering, for simplicity, the case of W production, we find that the symmetrization procedure 
of ps| ] simply amounts to replacing the n-dimensional 7" by the four-dimensional one. This clearly does not help 
remove the spurious UV terms; a counterterm will still be necessary. Of course, the ambiguity about whether to 
put the 75 on the right- or left-hand side of the 7" is no longer present. 

The string of gamma matrices ( |65|) gets contracted with a tensor I pa arising from the vertex loop integral. This 
is a good example to illustrate how the unphysical term arises in DRED. We will keep the argumentation general 
so that it is valid for both massless and massive quarks. Firstly, we have in DRED, 



I pa has the form 



y, 7 P 7 a 7 a^ = _27-7 a 7". (67) 



I pa =9 n pa C 2i + T pai (68) 



where T pa is a tensor which depends on the external momenta and C24 is a scalar coefficient. Performing the 
contraction gives 

= -2C 24 [(1 - eh a + 7?] + 7 CT 7 Q 7 P ?V, (69) 

where we used the relation, valid for n < 4, 

l°Tl p 9 P \ = -2(1 - eh a - 2 7 f . (70) 

It is precisely the term ~ 7" which is removed by the DRED counterterm. The remaining terms give the correct 
result, i.e. that required to satisfy the QED Ward identity between vertex and self-energy graphs. That identity is 
most straightforwardly checked by retaining the quark mass, in order to avoid mass singularities. Having satisfied 



(66), the virtual corrections are analogous to those of the 7-g vertex, since the Born term simply factors out. 

Using the Born term of the specific scheme, the net virtual contribution in a form valid for all consistent 
n-dimensional regularization schemes is (ab — qq, qq) 

/ ,~ab,F±B \ 

/ a<T mn,BB' \ ab,F±B / U /, s^-, a s „, n 

2 2tt 2 2 P, 



- — - 7+ — - 

e 2 3 e C F 



(71) 



plus possible extra soft terms which would cancel with opposite ones in the soft bremsstrahlung, as discussed 
above. Here we defined 



The result ([71]) follows directly from that given in 



VI. BREMSSTRAHLUNG CORRECTIONS 



The leptonic integrated 2^3 particle squared amplitudes have the form 

KA ab.F±B r 1 f2Tra s H 2£ \ .~.ab,F±B j-,ab,F±B /-on 

■M-rana^ = M 1 M 2 J M mn,2->2 F mn,2->V ( 73 ) 

where the factor cf is present only for ab = qq, qq. The M are given by 
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M qq,F±B M qg, r ^D _ KA9q,r 3=-D KAq 
Jvl mn,2^2 ~ Jvl mn.2^2 " Jvi mn,2^2 ~ Jvi mn,2^2 



<qg,F±B 



,gq,F±B 



iq q,F±B 



(74) 



and 



KAqq,F±B _ j^qg,F±B _ ^gq,F±B _ **qq,F±B f? ,x 
JVl mn,2—2 ~ JVl mn.2-*2 ~ JVl mn,2^2 ~ M mn,2-t2' \'°) 



For the qg subprocess in production, the r.h.s. of (|7^), ( f75| ) implicitly contain a sum over the various quark 
flavors into which the gluon may split to form the W . 
For the qq subprocess, we have 



pqq,F+B _ pqq,F+B _ S (1 + W ) 



Pi • kp2 ■ k 



(76) 



pqq 

mn 



F-B 



— pqq 

mn 



F-B 



s 2 (l+w 2 ) s(l + w) 
Px-kp 2 -k pi-q 



For the qg subprocess, we have 



(77) 



<-) , <-) 
p q g/g q ,f+b 



2w + 2 



P2/1 ■ k | ^ A n/m P£ g (w) 
s p 2 /i ■ k 



(78) 



(-) 

p q g.F-B _ 2 



pi ■ q 



2w + 2 



P2-k A n P^ g (w) 

h s — 

s p 2 ■ k 



(79) 



where 



(-) 
h q 9 

unfit 



-1 + w, 



(-) 
h q 9 

IL ll/ ul 



1 — W 



(80) 



and 



pg 

mn 



(-) 

q ,f-b 



= 2 



pi-q 



2w~2 



Pi-k A m P 4 H 

h S 

s p\ ■ k 



(81) 



with 



(-) 
q 



h 9 q , , = 2w 2 , 

uu/ui ' 



(-) 

Kn = 2w - 

11 1 Lu 



(82) 



Various checks on the above results are possible at this point. From crossing symmetry, it follows that 



M qg,F±B 



1 



2C F 



Mf u F±B ( P 2 - -*), Ml 9 u F±B 



1 



2C F 



Mt F±B iP2 ^ -k) 



(83) 



and 



MZ F+B = -±-M^ F + B {pi - -k), 



M gq,F+B 



2C, 



-MV U F+B { P1 - -k). 



(84) 



The above correspondences were explicitly verified. ( |84| ) holds only for the F + B case, since the F — B integration 
depends on pi and hence destroys the crossing symmetry. 

Also, the q <-> q interchanges were checked using CPT invariance. The only difference between the F — B and 
F + B cases is that in the F — B case, one picks up an extra minus under CPT, since CPT implies l + «-> l~ and 
the leptonic integration is antisymmetric. 

The 2 — > 3 particle bremsstrahlung phase space is given by 



da 



ab,F±B 
,BB' 



dM 



l-2e 



Br 



T(l-e) 



dyy e 0--yY 



VPz\M\i°?, BB , 



F±B 



(85) 



where 



1 + cos ( 



(86) 



and 9 is the angle between p\ and k in the p±, p2 cm. 



ft 



Define 



T ab,F±B 



_ (l-w) 1 - 2e r(l-2e) f 1 



Then, noting @, @, 



ab,F±B ' 
,.BB' 



\r„] — v ab ' F±B (F)in 1+e r(F)r ab < F±B 
27T " l "' S - B °l e / J mn 



-e pab,.F±B 



(87) 



where the Xmn ±B are defined analogously to the _M„ri F±B (see Q, (|75|)). The resulting integrations are rather 
straightforward. Using the approach of [Q we may cast the I'^^ ± m a form valid for all consistent n-dimensional 
regularization schemes, making use of the n-dimensional splitting functions A^PPj. As for the virtual graphs, we 
omit a possible extra soft piece which would cancel in the virtual plus bremsstrahlung sum and use the Born term 
appropriate to the regularization scheme in question. 
We obtain, for the qq subprocess, 



e 2 [ ' e C F 



ln(l — w) 
1 — w 



A(l + w)Hl-w) + k^ ±B (w), (89) 



with 



and 



■„qq,F+B 



(w) = -2(l-w), 



F—B 



(w) = 2(1 +tu) lnw 



A mn Pij — A m Pij -\- A n Pi^ 



For the qg subprocess, we obtain 

1 



„/,,f ±B = A ^p» W + 21n(l - w)A:!> /m P* g (w) + k a ± F±B (w), 



with 



l.q g,F+B 
(-) 

k i b,f-b 

uu : lu 
(-) 

h i g,F-B 

,gq,F-B 



h 9 



(-) 

q ,f-b 



1 1' 



(i - «0 

4 

(1 - «0 
4 

4 

(1 - w) 



(1 + 2>w) + (1 — w) In w, 
(1 + 3w) — (1 — w) lnw, 
(l + 3w)-2u; 2 Into, 
(1 + 3w) — 2iu hiw 



and 



a ab p 



A n P qg :b = g 
A m P qg : a = g 



(90) 
(91) 

(92) 

(93) 
(94) 
(95) 
(96) 
(97) 

(98) 



VII. FACTORIZATION OF MASS SINGULARITIES 

89] ) represents a soft divergence (and simultaneous mass singularity) and is cancelled by 



The term - 1/e 2 

an opposite term in (fn). The remaining terms ~ 1/e in (|7l]), (|89|), ( p2| ) represent mass singularities and do not 
cancel. Hence, they must be removed by expressing the bare parton distributions (and fragmentation functions 
in processes with final state mass singularities) in terms of the renormalized ones |30[ ]. Thus, our parton model 
expression ( |l8| ) is understood to be initially written in terms of the unrenormalized parton distributions, whose 
relation to the renormalized ones is given by 

c(e) x - atj ^dy rJ / A 



3 

x{A k P*(x/y) + eA k T l3 (x/y)}, (99) 



c ( t) \- uy / ay f r,j/A ( 2 ^ 
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^4(4.ri^4-^ +ln47r+0(£) (100) 



where 



and ctj — a s , unless j (=A) = 7, in which case ctj = a. After making the above substitution in (|18[), we recover 
the same form, except that now the subprocess cross section is finite and the parton distributions which enter are 
the renormalized (finite) ones. This will be demonstrated explicitly later. Note that photons, like hadrons, have 
internal partonic structure. Hence ( |99| ) is not restricted to hadronic initial states. We do not introduce the any 
additional mass scale here associated with the mass singularities as is commonly done. In this form, the above 
definition is consistent with the usual MS definition of renormalization when A^Ty = 0, as will be discussed 
shortly. 

The fragmentation functions T> A ^, which represent the probability for quark i to split into a collinear hadron 
(or photon) A, with momentum fraction z, have the following renormalization: 

v k ( Z 'V ) = v k ( z >v ) 



x[A k Pf i {z/y) + sA k T ji (z/y)). (101) 

All of the freedom in the factorization scheme is parametrized by the subtraction terms Tij[x/y). In theory, 
the Tij could depend explicitly on x as well. Although we will discuss constraints on the allowed Tjy in some 
detail later in this section, it is worthwhile to briefly summarize the most commonly used schemes here. In Table 
I we list various schemes of interest. For the MS (as opposed to MS) version of any of these schemes, simply add 
(je — ln47r)AfeP^ to the corresponding Afcly. 

The e modified minimal subtraction scheme (MS e ) was introduced in Q, the 'polarized' modified minimal 
subtraction scheme (MS P ) was introduced in |25| and the 'helicity conserving' MS scheme (MS/fc), as we shall 
refer to it, corresponds to the scheme used in determining the two-loop polarized splitting functions in pi]]— fl33| , 
although it was inaccurately referred to as MS there and in many other papers making use of it. This inaccuracy 
appears to have been acknowledged in j33] . 

In HVBM regularization, the MS scheme leads to helicity nonconservation of massless fermions. Hence it is 
not generally used. The other three schemes are helicity conserving, at least in the direct sense. A potential 
problem in the MS#c scheme will be discussed later in this section. For the MS#c an d MS P schemes, one 
must use the MS scheme as the corresponding unpolarized factorization scheme when comparing polarized parton 
distributions to unpolarized ones. This is implicit in the definition of those two schemes, which are defined only 
for polarized parton distributions. This is an advantage since unpolarized parton distributions determined in the 
MS scheme are widely available, in contrast to the MS e scheme. On the other hand, it is not clear what effect 
these two schemes could have on various positivity constraints which require polarized quantities to be smaller 
in magnitude than the corresponding unpolarized ones, as will be discussed briefly at the end of this section. 
Nonetheless, presently available NLO polarized parton distributions are predominantly determined in the MShc 
scheme, so we will use it in making physical predictions. 

The MS £ scheme is the only of the three helicity conserving schemes with the following properties, (a) It treats 
the polarized and unpolarized parton distributions analogously, (b) It gives regularization scheme independent 
analytical results and parton distributions (equivalent to those of DRED MS), (c) It satisfies the supersymmetric 
identity 

A k P gg + A k P gg = A k P~ gg + A k P gg , ,g = gluino, (102) 



for the unpolarized |3jj35[| and polarized [pl| , |32| two- loop spacelike splitting functions as well as the timelike 1 33 1 
splitting functions (both polarized and unpolarized). This is due to the equivalence with DRED, which in turn 
implies the app licab ility of four-dimensional calculational techniques. The timelike splitting functions are those 
which enter in ( |l0l[ ). They only begin to differ from the spacelike ones at the two-loop level. 

Substituting (|99|) in ( |l8| ) and convoluting with the n-dimensional Born term of the regularization scheme being 
used generates several factorization counterterms. Those associated with the qq subprocess (ab = qq, qq) have the 
form 

"^if^J = lx:t F BB^)^ 1+£ C(e) (±j [A mn P» + eA mn T qq (w)]. (103) 
The remaining terms are associated with the qg subprocess (ab = Q g,g q ) and have the form 

/ d(T mn -BR' \ I 



dM 



— ) = \xXBB^)^ 1+£ C{e) (£\ [A^P^H + eA% m T qg (w)}. (104) 

Ct 
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From Q we may infer the relation between parton distributions of any regularization and factorization scheme, 
denoted 1, with any other, denoted 2: 

= fil£(x) + E ? f-fffm^kPZaWv) - A k T hlJ (x/y)} 
j 2vr J x y ' 

-[AkP^jix/y) - A k T 2jij (x/y)]} + 0{a 3 a s ). (105) 
Similarly, for the fragmentation functions, 

^ 2tt J z y 

-[A k Pl ]t (z/y) - A k T 2di (z/y)}} + 0(a 3 a s ). (106) 

Going from scheme 1 to scheme 2 simply amounts to expressing the parton distributions and fragmentation 
functions of scheme 1 in terms of those of scheme 2, as was shown in 0. Also, we see that working in any 
consistent n-dimcnsional regularization scheme and any factorization scheme is equivalent to working in any 
other consistent n-dimensional regularization scheme with the appropriate factorization scheme. Thus, according 
to our definition of consistency given in section III, a consistent n-dimensional regularization scheme must (a) be 
mathematically consistent, for the class of subprocesses under consideration, and (b) have a process independent 
set of splitting functions within that class, in order that it may be related to the HVBM scheme as described in 
section III (assuming DREG and HVBM themselves do not fail for some process, unlikely as it may seem). The 
DRED scheme satisfies both criteria for one-loop QCD corrections to polarized processes. 

From j| we may also infer that the conversion term which must be added to the subprocess cross section 
calculated in scheme 1 to obtain the subprocess cross section in scheme 2 is obtained via 

Vconv = 22a ctt i{[A k P^(wi) +eA k T liij (wi)] -> e[A k T 2 ^{wi) - A k P^ tij (wi)] 
i 

- e[A fc T M ,M - AfcPf^M]}, (107) 

where wi is the argument of the splitting functions in the counterterm a c t,i- The conversion terms are generated by 
making the above substitution in all factorization counterterms for both initial and final state mass singularities 
in general. As well, one must express the couplings of scheme 1 in terms of those in scheme 2 for processes where 
coupling constant renormalization enters. Then, the above form is valid for all one-loop Q CD c orr ectio ns, not 
just Drell-Yan and DIS. Also, it is straightforward to verify that making t he s ubstitutions ( |105| ), (106) in any 



one- loop cross section calculated in scheme 1, leads to the conversion term (107). 

Parton distributions determined using different processes, but the same scheme, will differ by C(a™ +1 ) when 
C(a") corrected cross sections are used in the fits. Since a basic assumption of the parton model is that the 
all-orders parton distributions will be process independent, the process dependence of the parton distributions 
can be seen to arise from the neglected higher orders. The physical parton distributions are in fact the all- 
orders parton distributions since, in reality, it is the all-orders subprocess cross section that they are convoluted 
with. Hence, if the all-orders parton distributions were not process independent, there would exist no process 
independent parton distributions. By all orders, we mean sufficently high order in a s that higher orders may be 
neglected for all practical purposes. The process dependence is therefore easy to understand since some processes 
have large QCD corrections (i.e. Drell-Yan) while others have relatively small corrections (i.e. DIS). Hence the 
extent to which the 0(a") and C(a" +1 ) parton distributions differ, in any particular scheme, will depend on the 
process from which they are determined. So, for low n (i.e. n = 0), processes with significantly different QCD 
corrections will yield somewhat different parton distributions, but the ones with smaller QCD corrections will 
lead to parton distributions which are closer to the all-orders ones. One expects, on the other hand, that the 
ratio of the polarized to unpolarized parton distributions will not have a large process dependence since, even 
for Drell-Yan like processes, spin asymmetries tend to exhibit perturbative stability as we shall demonstrate in 
sections IX - XL 

Equation ( jl05| ) is very powerful and restrictive. This is because it has implications for various conservation 
rule motivated sum rules. According to the parton model, the number of valence (constituent) quarks, q v , in a 
hadron, A, is given by 

N*"' A = I dyf^ A (y) = const. (108) 
Jo 

More specifically, 

]\fd v /p _ f^d/p _ jyd/p _ 

N u v /p = N u/p _ N u/p = 2 ^ ( 1Q9 ) 



14 



The fact that this must hold for all energy scales, /x, is an expression of charge/probability conservation. Since 

N q l A is just the first moment of fu A {y) and since the k th moment of a convolution is the product of the fc th 
moments of the functions being convoluted, it is easy to see what happens to the above conserved quantity as 
one goes from one scheme to another. 

Firstly, we note that the gluonic contribution of ( |105| ) to the relation between quark distributions of different 
schemes cancels for the valence distributions. Suppose we wish to see what happens when going from rcgularization 



scheme 2 to regularization scheme 1. Using (105), we get, for valence quarks, 



K /A = N^{1 + g[(A„Pl, M - A u T hqq ) 

- (A u Pl m - A u T 2 . qq )]}, (110) 

where g denotes the first moment of g(y), for some function g. Suppose scheme 1 is DRED, for which 
^uPoRED.qq = 0, and scheme 2 is DREG, for which 

Au~P D REG m = 0, (HI) 

as can be seen from (p9|). Assuming the same factorization scheme, i.e. A u T2.ij = A u Ti t ij, we obtain 

N lJA _ N q v /A (u2) 
DRED ~ DREG' K 11 ^) 

so that the two regularization schemes are mutually physically consistent regarding charge conservation. The 
requirement, 

A u T qq = 0, (113) 

follows from the fact that the 0(a s ) corrections to charge conservation must be zero in n dimensions. By 
corrections, we mean that one can think of a zeroth order piece to A u P qq (z), A u P gg (z) which is proportional to 
5(1 — z) and denotes the no interaction scenario. Then, it is clear that the usual leading order splitting functions 
denote 0(a s ) corrections to the non-interacting case. 

Let us now consider momentum conservation. Let P l / A denote the total momentum fraction carried by parton 
i. It is just the second moment of fu A : 

P l ' A = fdy V rJ A {y). (114) 
Jo 

Of course, momentum conservation must be applied to the sum over all partons. One can easily verify that if, in 
some factorization scheme, the DRED parton distributions satisfy the momentum conservation rule 

J2P i/A = l, (115) 

i 

then the DREG parton distributions, in that same factorization scheme, will also satisfy it, and vice-versa. The 
basic reason is that the requirement, which DREG satisfies, that 

fdzz[A u P r qq (z) + A u P^ q {z)] 
Jo 

= f dzz[A u P gg (z) + 2N F A u P^ g (z)} = 0, (116) 
Jo 

(which follows from the fact that the 0(a s ) corrections to momentum conservation must be zero, in n dimensions) 
leads to cancellations in the sum over partons, so that the net difference in total momentum between schemes is 
zero. This is easily verified from Eq. (J29|). The individual quark and gluon momentum fractions will depend on 
the regularization (and factorization) scheme as well as the energy scale, fj,, to 0(a s ), but a change in the quark 
momentum fraction is cancelled by an exactly opposite change in the gluon momentum fraction. Thus, the total 
momentum is the same in both schemes, and is independent of /i, as follows from DGLAP evolution p6| , p3[ . 

All this discussion tells us that DRED and DREG are mutually physically consistent in general, in the un- 
polarized case. One can convert unpolarized parton densities from one regularization scheme to the other and 
the conserved quantities retain their values. The real problem arises when one converts from one factorization 
scheme to another. It is safe to go from MS to MS since, in the MS scheme, all the AfcTy are proportional to 
AkP^j (je — In 47i"), and hence do not change the conserved quantities. Also, since DRED is equivalent to the MS e 
scheme, we see that MS £ and MS are mutually physically consistent a s we ll, for unpo larized processes. In general 



though, for arbitrary A^Tij, there are no analogous constraints to (113) and ( |116| ), at least without invoking 
constraints from higher order evolution. This means that the aforementioned conserved quantities will differ in 
general to 0(a s ) when going from one arbitrarily chosen factorization scheme to another. Since it would be quite 
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a formidable task to measure so precisely the parton distributions that one could definitely verify all conserved 
quantities to high accuracy, how do we know which schemes are correct a-priori? This is even more pertinent to 
the case of the polarized parton distributions, where the data is definitely lacking. 

We may gain some insight into this issue by studying the scheme dependence of the polarized (nonsinglet) 
distributions. The net spin, S q / A , carried by quarks of flavor q, in hadron A, and the corresponding total quark 
spin, iAEyi, are given by 

sq/A = \ J/yft^' = E sq/A > 

respectively. The spin carried by the valence quarks is conserved under one-loop evolution due to chirality 
conservation. More formally, A;P gg = 0, substituted in the expression for the first moment of /; ', evolved using 
the DGLAP equations gives S q "/ A — const, under one-loop evolution, as was the case for N q "/ A . For a discussion 
of two-loop valence evolution effects, which in general violate the above conservation, the reader is referred to $7\ . 
The main point of interest concerning two-loop evolution is that chirality conservation is relevant to nonsinglet 
combinations of the form ff + f q — f? — f q , and linear combinations thereof. The two interesting ones A^3, 
AA$ satisfy 



where 



AA 3iS = const, (118) 



AA 3 = f? + ft - ft - ft 

A^s = ft + ft + ft + ft ~ 2(/f + ft)- (119) 
Similar argumentation to that used in relating N{° to N% leads to 

A% fi = AAl fi {l + ^[(&iP E i, qq ~ A,Ti, M ) 

- {^iPl m - A ; T 2 , W )]}, (120) 

From Table | and Eq. (|o|) we get 



—DREDQAS) _ -jHV(MS E ) _ —HV(MS HC ) 
^^3,8 — ^^3,8 — ^^3,8 

Thus, the MS scheme of DREG(HVBM) gives different values for A^3, AAg than all the 'chirality conserving' 
schemes. This is an explicit example of the factorization scheme ambiguity. 

The real implication of this seeming nonuniqucncss problem is that certain factorization schemes may be ruled 
out from the beginning since they arc manifestly unphysical. One can determine which schemes are unphysi cal b y 



imposing the fact that the unrenormalized parton distributions ( p9\j themselves satisfy the conservation rules ( 108), 



(115) and (118) and yield the correct value for the corresponding conserved quantities (except as will be mentioned 
below). The reason for this is that if a particular regularization preserves the necessary symmetries: unitarity, 
translation invariance and chiral symmetry, then the unrenormalized parton distributions in that regularization 
will obey the corresponding conservation rules: conservation of charge/probability, momentum, and massless 
quark chirality. Hence, the connection between renormalizcd and unrenormalized parton distributions must 
not destroy those symmetries in the renormalizcd parton distributions. If, on the other hand, a regularization 
violates a particular symmetry, a finite AfcTij will have to be introduced to insure that the renormalized parton 
distributions obey the corresponding conservation rule, since the unrenormalized parton distributions themselves 
will manifestly violate it. 

Thus, using similar rationale to that used in going from one regularization scheme to another, we arrive at the 
physical consistency constraints on allowable factorization schemes, 

dzA u T qq (z) = 0, (122) 

/ dzz[A u T qq (z) + A u T gq (z)] 
Jo 

= [ dzz[A u T gg {z) + 2N F A u T qg {z)} = (123) 
Jo 

and 



1G 



^T qq (z) - A,i* (z) = A u T gg (z) - A u PUz) 



(124) 



The reason we had to subtract A^P^ in (124) is that for regularization schemes like HVBM which violate chi- 



ral symmetry (AiP qq ^ A u P qq ) 7 the unrenormalized parton distributions themselves manifestly violate helicity 
conservation. One must always subtract the helicity nonconserving part in AiP qq in order to restore this conser- 
vation. This is why HVBM(MS) disagrees with the other schemes and is manifestly unphysi cal. Under two-loop 
evolution, helicity non-conservation in HVBM (MS) manifests itself via direct violation of ( |11§| ). This may be 
seen by taking scheme 1 to be HVBM(MS) and scheme 2 to be one of the 'chirality conserving' schemes, then 

differentiating (|l20|) with respect to In /i 2 . Comparing with the DGLAP equation for A A 3 8 one finds precisely 



the same nonvanishing result for dAA 3 8 /<iln/i as follows from the two- loop analysis 31 32j|. 
The condition (124) is somewhat stronger than the conservation rule (118). It also insures that 



<r(helicity flip) = 0, 



(125) 



where cr (helicity flip) generically denotes any doubly polarized (in the initial state) subprocess cross section where 
a massless quark flips helicity, with the helicity at both ends of the quark line fixed. 

The relat ion between the total quark spin in scheme 1 and that in scheme 2 is obtained by taking the first 
moment of ( 105 ) and summing over quarks (including antiquarks), 



2N f f°[(A l P[ tqg - A,T hqg ) - (AiP 2tQg - A ; T 2i9ff )] 



(126) 



Since DRED preserves chiral symmetry, it is natural to use it as a reference scheme. If we also desire that the 
total spin of the quarks be the same as in DRED(MS), then AiP qg ^ in HVBM regularization leads to violation 

of this in the MS^c scheme. We must impose 



I dz[AiT qg (z) - AjPg (z)] = 
Jo 



(127) 



to restore equivalence. Hence, of the factorization schemes listed in Table Q, only MS e and MS P give the value of 
total quark spin equal to that in DRED (MS). 

One may ask if there is some deeper physical reason why the MSijc scheme gives a different value for the total 
quark spin. We notice, using (||) and pi), that in the HVBM scheme 



qg, 



(z)^P"_ +(l-z). 



(128) 



This means that when a polarized gluon splits into a collinear qq pair, the q and q do not necessarily have opposite 
chiralities, as required by quark chirality conservation. The above result is possible since, if there is an initial 
state gluon, (125) does not apply. 

The MS e and MS P schemes correct ( |128| ), while MSirc does not. More precisely, strictly speaking, chirality 
conservation imposes the indirect constraint, 



[A u T qg (z) + AiT qg {z)\ - [A u i* (*) + A,l* (*)] 

= [A u T qg {\ - z) - A,T 99 (1 - *)] - [A„P ? e g (l -z)- Ai?;(l 

=> A t T qg (z) - AtP^z) = -[A,T m (l - z) - A l Pl g {l - z% 



(129) 



which both MS £ and MS P satisfy, but MS^c does not. The reason, as usual, is that it is the difference, A^P? (z) — 
AfeTjj(z), which enters in physical cross sections. Thus, th e su btraction term, A^T i j(z), effectively modifies 



the splitting function, AkPh{z). The unpolarized parts in fll29| ) cancelled due to the unpolarized constraint 
A u T qg (z) — A u P qg ( z) = A u T qg (l — z) — A u P q (1 — z), which is satisfied in all the schemes considered here. 



Satisfaction of (129) implies satisfaction of (127), while the r evers e is not true. An interesting analysis which 
arrived, from a different perspective, at a conclusion similar to (|l27| ), in the context of DIS, may be found in [j38f . 
How strictly one wishes to impose chirality conservation is perhaps a matter of taste, however. 

Under one-loop evolution, the total quark spin is conserved. Under two-loop evolution, conservation or non- 
conservation depends on the choice of factorization scheme. Analogously to the nonsinglet case, we see that if AS 2 
is conserved and ASi ^ AE 2 , then AEi w ill n ot be conserved under two-loop evolution. The scale dependence 
of AEi may be obtained by differentiating (126) directly. The MS £ and MS P schemes lead to conservation of AS, 
therefore MS#c does not. Alternatively worded, AS is not conserved in the MS^c scheme and the difference in 
the schemes exactly compensates this non-conservation such that AS is conserved in the MS e and MS p schemes. 
The evolution of AS in the MS#c scheme obtained via differentiation of (126) agrees with the two- loop analysis 
|3l|j32]| and can be seen to arise from the non- vanishing of AiP qg . Under two- loop evolution, it is the non- vanishing 
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of the first moment of AiP qq (from the pure singlet part) which leads to the non-conservation of AS, however. 
Hence the non- vanishing at two- loops of AiP qq is seen to be a consequence of the non- vanishing at one- loop of 
AiP qg , while AiP qg remains zero at two l oops , in all three schemes. 



At any rate, we see that satisfaction of (127) leads to conservation of total quark spin. The DGLAP evolution 
equations will be elaborated on in the next section in connection with scale dependences. It is worth noting that 
analogous argumentation can be used to relate the spin carried by the gluons in different schemes. Since such 
differences do not enter in the observables considered here, at one loop, we will not discuss them. Also, direct 



differentiation of (105) gives the correct scheme transformations for all two-loop splitting functions, generalized 
to include the regularization scheme dependence (i.e. Ty — > Tjj — iy). 

So far, we have not discussed the interesting requirement of positivity of the parton distributions, 

[i 2 ) > 0, fjj + (x, M 2 ) > 0. (130) 

The only guarantee of positivity applies to the physical cross section itself, as measured by experiment. Individual 
subprocess cross sections need not satisfy positivity due to the interconnection of the subprocesses arising from 
the renormalization of the parton distributions discussed earlier in this section. To the extent that leading order 
cross sections dominate and that the parton distributions are process independent at leading order, positivity 
should not be a problem. On the other hand, for processe s wi th large radiative corrections, some factorization 
schemes could lead to parton distributions in violation of ( |l30| ). This issue has not been thoroughly studied at 
this point, but it will become important as more experiments involving polarized hadrons are undertaken. 



VIII. ANALYTICAL RESULTS TO ONE-LOOP ORDER 



We may now present the analytical results for da^^^g, / dM in a form valid for all consistent n-dimensional 
regularization schem es an d for a general factorization scheme, starting with the qq subprocess. Combining (j57|), 
(71), fl8S|), (|89| ) and (103), we obtain the total result (ab — qq,qq), 



da, 



ab,F±B ' 
mn.BB' 



dM 



NLO 



xtn F BB> (0) *(* ~ «0 + \(\ - 7 ) ' ' ' - " ' 



ln(l — w) 
1 — w 



AmnP^jw) 

C F 



In- 



M 2 

U>/1 2 



4(1 + w) ln(l - to) 



+ — [A mn T qq (w) - A mn pe g (w)) + k<£ n F±B ( 



C F 




where the fc„^ F±B are given in (|90|). 

We now present the result for the qg subprocess. There is no 0(1) term (in a s ). Combining 

(-) (-X 



(131) 
I), @ and 



(104), we obtain the total result (ab = q g,g q ), 



' ,~ab,F±B 
aa mn,BB' 

dM 



~ab,F±B fr,\ a s I » a b r>4 / \ 



NLO 



M 2 

In ? + 21n(l - w) 

WfJ, 



[At, m T qg (w)-Al) m P^ qg (w)\ 



„ab,F±B 



(w) 



(132) 




where the k^^f ±B (w) are given by (|92 

We see at this point that the results of W (for the uu and //, F + B. vector coupling case with massless 
B = B') are contained within the results presented here. Those, in turn, were found to be consistent with various 
existing results calculated in certain specific schemes, for both the uu ||{],|l| and 11 0J^] cases. Single-spin W ± 
production in hadron-hadron collisions has been studied in Considering the subprocess ab — > W ± c, where a 

is longitudinally polarized and b is unpolarized, we see that for ab = qq, qq and for a = g,b = q , our expressions 

reproduce those given in Appendix B of ||. For a = q ,b — g, it seems at first sight that there is a discrepancy. 
In order to reproduce the result of ||, one needs A u P^ g = —1/2, in disagreement with (|2^). The reason for 
this disagreement is that Q averages over 2, rather than 2 — 2s, gluon polarization states, as can be inferred 
from their Eq. (10). That convention does indeed lead to A u P qg = —1/2. Unfortunately, all of the presently 
available unpolarized parton distributions are determined using the convention of 2 — 2e gluon polarization states. 
Consequently, the expressions ofj3[, for the qg subprocess, should be appropriately converted using the technique 
of [|]], which amounts to Eq. ( |107| ) of this paper, before being used in numerical studies demanding complete 
consistency at NLO. 



18 



At this point, we may investigate the scale dependence of the cross section. Define 



da 



AB.F±B 



dM 



(133) 



and let a denote the contribution to a which depends on [i 2 at 0(a s ) coming from a — q, b — q; a = q, 
b = g and a — g, b — q for one quark flavor and for production by one bosonic interference channel. These 
subprocesses are connected via the renormalization of the parton distributions, but are disconnected from the 
remaining subprocesses in that sense. We have (taking = = 1) 



or = y^ F±S 
Aran 



(0) 



dx„ 



dXb,afm( X a, ^)fn( x b, f^)5(l ', 

T / X a,b %a%b 



-Hi f dx b -mx b ,» 2 ) C ^[/«(w 2 )A m P 9 4 >)ln/x 2 + /s( W 2 )A m P 9 4 g Hln M 2 

Z7T J T Xb J r/x b X a 

~ C dx a -fl(x a ,^) [ ^[/^^^A^V^^^ + Z^^^^An^H^M 2 ] 



/x a X b 



(134) 



where 



fk( x a,b, M 2 ) = fk ' ( x a,b,H 2 )- 



(135) 



In DIS, one has an analogous expression, except without the convolution with the term. The independence of 
a on /i to 0(a s ) implies 



da 



d In /i 2 



Here we note 



da s 
din fi 2 



47T 



0{al). 

a 2 s = 0(a 2 s ). 



(136) 



(137) 



Taking into account the fact that the evolution of fm A is independent of fn B and fn , we may arbit raril y 
set, for any fixed fi 2 , /^(xh,/i 2 ) = and f%(x b ,n 2 ) = (5(1 — x b ), after performing the differentiation in (136). 
Of course, fn B and fn B still evolve with /i 2 . We then use (136) with m = n, / r f(o; a ,/i 2 ) = /^(xf,,^ 2 ) = 
and /^(x a ,^ 2 ) = (5(1 — x a ), f^{x bl ^ 2 ) = (5(1 — x b ) to solve for and eliminate df%(x b , fi 2 )/d\nfi 2 . The result for 
dfl{x bl ^i 2 )/d\n^ 2 is simply ^A n P* q (x b ). 
We thus obtain 



d fk( x ^ 2 ) = Usjy 2 ) f 1 dy 
dln/i 2 2tt 



■ir k (y,H 2 )A k P^(x/y) + f k {y, ^)A k PUx/y)] + 0(a 2 s ), 



(138) 



for k = u,l. This is the DGLAP equation |36J,|23J with one-loop splitting functions. Therefore, one needs at least 
one-loop evolution to guarantee \j? independence to 0(a s ). Gluonic evolution only enters at the 0(a 2 ) level. 
We may write, using the Taylor expansion for fixed x, 



tf(x, f i 2 ) = f q k (x, t i 2 ) + 

= f q k (x,f4) + 



2tt Uo 



2tt 



-i 2 



C i,fc( a '> A*o) 
c i,fc( a; ' A*o) 



/«(a;,M 2 ,/i 2 ), 



2tt 



In 



c l,fc( x 'Mo) + 



(139) 



where 



c l,k( x ' A*o) 



(140) 



One can write an analogous expression for /|(x,/i 2 ) using the DGLAP evolution equation for f k (x,/j, 2 ), yielding 



Jx y 



(141) 
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Then one could express a in terms of the / and parametrize explicitly the 0{a 2 s ) fj, 2 dependence with respect 
to some fixed scale /Zg. Suppose /Xq is a typical scale of the process (i.e. M 2 ), then there will not be any large 
logarithms and ct(/Xq) will have reasonable 0(a s ) corrections when /i 2 = /ij > A 2 . Having identified /j,q as being 
a scale which leads to reasonable 0{a s ) corrections, we now ask: what happens when we vary fj, 2 away from /x 2 ? 
Using (13E) in the parton model expression, we see that there will be /x 2 dependent terms will be proportional to 



Zir 



ln(^ 2 /A 2 ) 



(142) 



times an overall factor having roughly the same magnitude as the Born term. For typical values of \i and a s , 
ln(/i 2 / Ho) is a small number so that the /j, 2 dependence is rather suppressed so long as /s 2 is not too far from 
Hq. These terms will only cancel when two- loop corrections to the subprocess cross section are included. There 
are also terms with lower powers of ln(/i 2 //!Q) (i.e. from (137)) which require two-loop corrections for cancellation 
as well. 

The next question is: how much better are we doing than a Born level calculation? In t he Bor n level result, 
the /j, 2 dependence comes entirely from the evolution of the distributions. Hence, from (134), ( 139 ), 



C Born — CBorn 



2tt 



(143) 



again times an overall factor having roughly the same magnitude as the Born term. We see explicitly that 
inclusion of one-loop corrections reduces the /i 2 dependence by an extra factor of ^ \u(/j, 2 /iiq) for reasonable 
values of ^ 2 . The above arguments apply to all one loop calculations in hadronic processes. In general, though, 
there is also a contribution from the evolution of the gluon distribution at 0(a s ), relative to the Born term. Thi s 
is because, in general, the gluon distribution may enter in the Born term. Then, imposing the analogue of ( |l3(j ) 
would determine c\ k as given by (141). 

One could also convolute the Born cross sections with parton distributions having no evolution (i.e. determined 
at some fixed scale jl 2 ), which would correspond to an "off-shell" renormalization of the parton distributions in 
higher order. But then the physical predictions would be wrong by terms ~ #^ p 2 ) 1 since the one loop 



corrections would contain such terms (analogous to the terms ~ ln(M //i ) present in ( |13l[ ), fll32j )). So, even 
though there is no renormalization scale dependence, the error will be larger for fio ^ P- or Mo "C ft- One expects 
an analogous tradeoff when using one-loop evolution instead of two-loop evolution in NLO calculations. 

Actually, this argumentation is somewhat naiive in the sense that we only considered explicitly the ~ 
[w= ln(/i 2 //i§)] 2 behavior. For processes such as the Drell-Yan process, which receive large tt 2 corrections, there 
will be residual /^-dependent effects arising from the ;U -dependent part of the parton distributions multiplied by 
the large correction terms in the cross section which will not cancel at the one-loop level. Thus, there will be an 
additional fj, 2 dependence of the cross section proportional to 



2tt' 

7T 



2tt 



ln(/z 2 //x 2 ) 



2tt 



(144) 
(145) 



which is still suppressed relative to the /i 2 dependence of the Born term for typical energy scales. Of course, there 
may be other sources of large corrections not strictly proportional to tt 2 (see the discussion near the beginning 
of section XI), so our tt 2 factor is really a generic representation of large correction terms in general. The rule 
then is that the larger the one-loop corrections, the greater the fj, 2 dependence. Inclusion of two-loop corrections 
(and evolution) will reduce this factor. We will not discuss two-loop effects in any detail here. Also, we will not 
perform detailed numerical studies of scale dependence, as this has been done extensively in the literature. The 
above discussion was included simply to clarify the origin of the scale dependence and to show how it is reduced 
by the inclusion of QCD corrections in general. 



IX. DRELL-YAN IN LOW ENERGY RUNNING AT RHIC 



Throughout, for unpolarized cross sections, we use subprocess cross sections determined in the MS scheme of 
DREG. They are convoluted with the unpolarized MRSG set |4lJ] of parton distributions. Varying the choice of 
unpolarized parton distributions produces negligibly small changes in the spin dependent asymmetries in general. 
Hence, we only consider the effect of varying the polarized parton distributions even though, strictly speaking, 
one should use specific unpolarized sets with the various polarized sets. We use the GSA and GSC sets p2[ as 
well as the GRSV (standard NLO MS) set E^] and compare the corresponding predicted asymmetries. Since all 
these sets are determined in the MS#c scheme, using HVBM regularization, we use subprocess cross sections 
determined in that scheme for doubly polarized predictions. Clearly then, for single-spin cross sections, we use 
the MShc(MS) scheme for the renormalization of the (un)polarized parton distributions which enter (i.e. for the 
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Tij). Consequently, the A fe P?- are always those of DREG/HVBM, given in Eq's (|29|) and (|30|). Throughout, we 
use a two-loop evolved a s , with 5 flavors. We also use /i = M in all numerical calculations. 

We will always present both the LO and NLO results for the asymmetries given. We use NLO parton distri- 
butions, polarized and unpolarized, in all asymmetries however. The reason, as discussed in section VII, is that 
the LO parton distributions possess a process dependence of order a s , while the NLO ones are process dependent 
only at order a 2 s . Therefore, it is in general not meaningful to use LO parton distributions in making predictions 
for processes with different structure from the one where they were determined; in the polarized case, DIS. In 
this way, our LO and NLO asymmetries are closer to the actual first two terms in the "all-orders" expansion of 
the asymmetry in a s and the difference in the LO and NLO predictions meaningfully measures the effect of the 
HOC. Finally, if one uses the LO polarized distributions determined from DIS, one would also have to use LO 
unpolarized distributions determined entirely from DIS as well, for consistency. 

In order to understand the qualitative features of the numerical results, in pp collisions, it is instructive to 
examine which values of x are being probed in the parton distributions for a given M. We will consider the 
leading order contributions for simplicity. From (p"8|), we see that da/dM involves an integration over xp (= xfq), 
equivalent to a single integration over x a (or Xb), both constrained by (|l^) with w = 1. da/dMdxp can be obtained 
from (^3|), though, and involves no integration. Numerically, it is found to generally peak at xf — 0. This implies 
that dominant contributions come from the region x a ~ Xb ~ \pr. The peak is not very sharp though, so this 
is somewhat of an approximation. Nonetheless, this feature holds for the energies considered here and becomes 
more pronounced with decreasing r. The basic reason is that at x a = Xb — yfr there is an equal contribution 
coming from the qq and qq subprocesses. As one moves away from xp = 0, one of the contributions becomes 
smaller since the q distribution is being evaluated at larger x. This is somewhat compensated by the other term 
getting larger, especially at larger r. Still, the net effect is a smaller contribution as if moves away from 0. As 
we go the lower t, the fact that the peak is more pronounced is not so surprising due to the 1/ (x a + Xb) factor in 
da/dMdxp, which is not present in the rapidity differential cross section. The rapidity differential cross sections 
generally peak at or near zero rapidity (x a — Xb) as well, although the peak may be rather broad, xp is more 
physical though, since it is simply x a — Xb at leading order (i.e. it is linear in x a and Xb unlike the rapidity). Hence, 
in what follows, including after QCD corrections, we will assume that the qualitative features of the numerical 
results can be roughly described by taking x a — Xb — \fr. Of course, all figures are obtained using the exact 
expressions only, since the approximation does not give quantitatively accurate predictions. We will also use the 
shorthand notation, 



with x ~ yfr. Also, 



,g^ft 9/p (x^ 2 ), Aq^Ag^ft'^ix,^), (146) 



da%f+ B /dM, (147) 



where a^ B was defined in (|l2|) as being the usual leptonic integrated cross section. 
Consider the double-spin asymmetry 



A " s vk*7t: (148) 

With the above mentioned assumptions, An is roughly given by 

(4/9)AuAu + (1/9) Ad Ad + {1/9) As As 
11 ~ (4/9)uu+(l/9)dd+(l/9)ss ■ ( ' 

Due to the extra factor of 4 and the relative largeness of the u v and Au v distributions, the cross section and 
asymmetry are mostly dominated by the up quarks. 

The polarized valence distributions are reasonably well constrained over a large range of x from polarized 
DIS. The polarized sea quark distributions, however, are only moderately constrained at smaller x, where their 
contribution in DIS is non-negligible. At moderate and large x, there is no constraint experimentally on the 
polarized sea quark distributions, as they are vanishingly small compared to the valence ones. 

Most models, however, assume the behaviour 

\Aq v \/q v X ^ 1, \Aq v \/q v 0, (150) 

which appears to be consistent with the DIS data. The result is that the asymmetries, and the distinction between 
the various polarized sets, will be largest at large r (i.e. large x) where the cross sections are small. Conversely, 
where the cross sections are large (small r) we expect small asymmetries if the polarized sea quark distributions 
are small at low all the sets assume. 

Fig.|l](a) gives the (virtual photon dominated) unpolarized Drell-Yan cross section, da/dM, at \/~S = 100 GeV, 
versus M for pp collisions relevant to RHIC. Here, we include only one type of lepton pair (i.e. ) . Throughout, 

we present the leading order, next-to-leading order with qq contributions only, and full ncxt-to-leading order 
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predictions for all cross sections and asymmetries. This means that, for the asymmetries, the numerators and 
denominators are treated in the same way with regard to which corrections are included. Here and throughout, 
we note that the qq corrections are positive and large. The qg subprocess makes a small negative contribution 
though. This highlights the fact that one cannot think of the qg subprocess as being physically separate from the 
qq subprocess. They are both related via the renormalization of the quark distributions. We will now investigate 
this issue thoroughly. 

The negativity of the qg contribution also holds in the MS e scheme. This will (slightly) simplify the process of 
understa nding the origin of the negative qg contribution. The analytical result for the subprocess cross section is 



given in ( 132 ). Since, when using MS e factorization, the scheme dependent part of the subprocess cross section 



is zero and since k'j^^ P+B (w) > and since fj, 2 = M 2 here, the only negative contribution comes from the term 
~ A u P£ g (w)ln(l - w). From @, (§?]), (§|), @ we see that this term arose from (1 - w)~ 2e A u P* g (w)/e and 
is therefore of collinear origin. This means that such terms could be factorized via the renormalization of the 
quark distributions, but would most likely lead to unphysical distributions violating some of the conservation 
rules mentioned in Sec. VII. So, we see explicitly that the issue of the negativity of the qg subprocess is intimately 
related to the issue of the renormalization of the quark distribution, which in turn connects the qq and qg 
subprocesses. 

The fact that the net qg contribution is negative tells us something about which regions of the w integration 



are giving dominant contributions. In order that the qg contribution be negative, we see from (132) that we must 
satisfy 

A u P g 4 g (u;)[21n(l -w) -hiw] + ^~ W \ l + 3w) < (151) 

(again working in the MS e scheme for simplicity). This is satisfied for 

w > .572022 . (152) 

Since the w int egration includes w = 1, we expect large contributions from the term ~ ln(l — w) from w near 1. 
In light of ( |l52j ), we can see clearly why the corrections are negative. 

The expected yearly integrated RHIC pp luminosity at V~S = 100 GeV is C = 160 pb~ (assuming it is linear 
in y/S down to this energy). This means that the statistics will be quite poor beyond M = 20 GeV, concerning 
asymmetry measurements, whose error goes like (see (pQ)) 



l-Al 



mn ^, A — Umn 



ly ^m^n V V & m C> n a uu O m O n 

where the approximate equality holds for A mn not too large. Here V is the number of events. Assuming 1 GeV 
binning, this means that for M ~ 5 GeV, A An ~ A%/T A T B . Unfortunately, this is a rather low mass scale 
where the parton model may not work well. For M ~ 10 GeV, AA U ~ 1.&%/V A V B . At low M (< 10 GeV), one 
also has to be careful of resonance backgrounds. Experimental cuts will make the errors slightly larger. On the 
other hand, with two independent experiments (at the PHENIX and STAR detectors), the combined asymmetry 
measurements should have roughly the errors given here. A more detailed error analysis is beyond the scope of 
this paper. 

Fig. |l|(b) gives the corresponding An. The behaviour is exactly as expected from the previous arguments. It the 
low mass region, where the statistics are good, An is of the order 1-2%. There is also little distinction between 
the various sets. This is a result of the common assumptions 

\Au\/u < 1, AS < (at small x). (154) 

There is certainly more freedom in Aq than is manifest among the various sets. The DIS constraints on Aq are 
rather weak and it will be interesting and important to see whether or not the low mass Drell-Yan asymmetry is 
as small (and positive) as it is predicted to be. 

It is also of interest to study the effect of the QCD corrections on An- From Fig. |l|(b), we observe good stability 
in the qq contribution in general. This is not unexpected, taking into account helicity conservation. For the GRSV 
and GSA sets, Ag is positive and sizable. For the GSC set, Ag is positive at small x, but negative and small 
at large x. Hence, the gluonic contribution is small and uninteresting in general for the GSC set. Therefore, we 
study throughout the gluonic contributions (and features in general) only for the GRSV and GSA sets. We notice 
that the gluonic contribution to An is always positive. The gluonic contributions to au and a uu are dominated 
by the ug subprocess. Since both Au and Ag are positive, the essential difference between the sign of the gluonic 
contribution to an and a uu is that arising from the sign of the respective subprocess cross sections. The most 
important difference between aff and is the difference in the overall factors (which are proportional to the 
respective Born terms, see (|73|)-(f75|)). There is a relative minus in the overall fact ors for and a u 9 



More precisely, for the correction, a u 9 , to be positive, we must satisfy (151) with A u P^ g (w) replaced by 
AiPq g (w). It is satisfied for 

w < .208233, w > .652396 . (155) 
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Using the same logic as before, dominant contributions will come from near w = 1, where the corrections are 
positive. Hence the corrections to an will be positive. The qg corrections to an are relatively larger then those 
to a uu since, for typical x and /i 2 , 

\Ag/Aq\ > g/q, (156) 

in the sets considered. Hence the qg subprocess typically weighs in more heavily versus the qq subprocess in an 
than in a uu . Thus we can clearly understand the sizable and positive qg corrections to An. This rationale also 
applies to the corrections to An in Z- and W /± -boson production. 

Figures |^(a) and ||(b) show the corresponding cross sections and asymmetries at \/~S — 200 GeV. All the 
features are the same. The luminosity is doubled and the unpolarized cross sections are larger for the same 
M (due to the smaller x, where the sea quark distributions are large). The asymmetries are somewhat smaller 
though. Hence the statistical significance (AAn/Au) is comparable. The main difference is that we probe lower x 
than at \/S — 100 GeV. Running at more than one energy is definitely an advantage in that one covers a larger 
range of x and /x 2 , compared with only running at one energy. This also allows some degree of cross checking 
via perturbative evolution and thus allows detection of various types of systematic errors. The net result is lower 
overall errors on the polarized parton distributions so determined. 



X. Z-BOSON PRODUCTION AT RHIC 



Z-boson production at RHIC is quite useful since it allows us to measure the polarized parton distributions 
at relatively large x. Also, the parton distributions enter as a different linear combination in the Z double-spin 
asymmetries than for production by virtual photons. This is helpful for disentangling the various contributions. 
One can also make use of parity violation to consider additional asymmetries which vanish for virtual photon 
production. The major limitation is the event rate. Certain asymmetries will be able to somewhat overcome this 
limitation, however. 

Fig. |^(a) shows the unpolarized cross section, da/dM for 5 < M < 125 GeV at \/~S = 500 GeV. Again, we 
consider production of only one type of lepton pair (i.e. The QCD corrections behave anaolgously to 

the virtual photon case. Using the expected yearly integrated luminosity of 800 pb" 1 and integrating the cross 
section between 80 < M < 100 GeV gives approximately 8,000 pairs. Hence we can measure asymmetries 

with an uncertainty (see (p0|), ( |153[ )) 

J- 1.5% 

AA mn ~ . B , (15') 

depending on the experimental cuts. Higher energy running would increase the Z event rate, but somewhat lower 
the asymmetries, as will become clear. 

The double-spin asymmetry in the Z-pole region goes rougly like, 

Au — ^ , 2 — ; — 2 x - : (158) 

l^qigiq + 9v q m 

evaluated near x ~ y/r = .18 . Putting in the appropriate couplings gives 

_ .29AuAw + .37(AdAJ + AsAs) 
11 ~ .29uu + .37{dd + si) ' ^ ' 

In general, for Z production, the GSA and GRSV sets correspond to two extreme solutions, with the GSC set 
lying somewhere in between. So, as in the last section, we will try to understand qualitatively only the GSA and 
GRSV predictions, since one can make a clear prediction for those sets. 

All sets assume an SU(3) symmetric polarized sea and take a negative down valence and a positive up valence 
distribution. For the GRSV set, the sea quark distributions are negative everywhere, while for the GSA set they 
are positive at intermediate and large x and become negative only at small x. Thus, for the GSA set, Aq > in 
the x range of interest, x ~ .18 . 

There will be some cancellation between the up- and down-quark contributions, but since 

the u contribution will still be bigger. Due to the smallness of the polarized sea quark distributions, Aq ~ Aq v . 
Hence, we expect 

An > : GRSV, An < : GSA (161) 

in the Z-~po\e region. 
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Fig. H(b) presents An for 5 < M < 125 GeV. The higher order corrections (HOC) have the sa me b ehaviour as 



for virtual photon production, for exactly the same reasons. As well, the sign of An predicted in (161) is verified. 
There is roughly a 4% variation in An between the two extreme cases. This means that one can rule out one case 
or the other, but not much more. Hence we must examine other asymmetries in order to do better. It would also 
be interesting to look at An in the low mass region, where the virtual photons dominate, due to the large event 
rate. As mentioned in the previous section, An may turn out to be larger than expected in that region. 
Define the single-spin asymmetry as 

A^A ul = ^=A lu = ^. (162) 

' u uu i u uu 

There is an overall minus relative to the definition often used [Q, owing to ([lj]), (|l5|). This asymmetry is nonzero 
due to the parity violating Z vertices. It is approximately given by 



.4, 



J2 q { 2 9aq9vq(q^q ~ ?Ag)] 

+ 9 2 vq ) (299)] 

Am + g a d9vdd v Ad 

(,9au + 9vu) UU + {9 2 ad + 9ld)( M + ss ) 

gau9vuuAu v + g ad gvddAd v 

OL + 9L) u u + (gL + 9ld)( dd ~+ s5 ) 

.lu v Au + .17d v Ad AuAu v + .17dAd v 



.29uu+ .37 (dd + ss) . 29uu + .37(dd + ss) ' 

where we took (A)q(x) — (A)q(x) as is done in all the sets considered. Noting that 

d(A8) 
u(.18) ~ 



(163) 



(164) 



and taking into accou nt (|l6C| ), we observe a rather large cancellation between the uAu v and dAd v contributions 



in the second term of (163). The net effect is that the first and second terms become comparable in magnitude. 
Hence Ai is directly sensitive to Au and Ad. Also, the Au and Ad contributions are comparable, unlike in 
virtual photon production. The second term is reasonably well constrained from DIS, while the first term is 
almost completely unconstrained (except in maximum magnitude). From the respective signs of the sea quark 
distributions, we expect Ai to be most positive for the GSA set and least so for the GRSV set. 

Fig. ^(a) shows A\ for 5 < M < 125 GeV. We observe the predicted behaviour. At small M, A\ vanishes as 
expected since the parity violating Z contribution also vanishes. The peak is just above the Z-pole and there 
is a well measurable separation between the various sets. So, the sensitivity to the sea quark distributions has 
improved as compared to An- 

One can take advantage of having two polarized beams to improve the statistical significance of the parity 
violating asymmetry. Define the two-spin parity violating asymmetry as, 



a pv _ cr (+,+) - «■(->-) 



^ 2(a ul /V B + a lu /V A ) 
p a,b =1 2(a uu + <j u /V A V B ) [ ' 



<r(+, +)+*(-,-) 

for V A ~ V B = V, (166) 



2 °»' r.„. '-n nH 



where the last approximate equality holds since an <C o uu as an involves two polarized parton distributions and 
is thus relatively rather suppressed. We see explicitly that A^ v is proportional to 1/V (relative to experiment), 
not 1/V 2 as is often assumed. In order to get an idea of the statistical error on Af t v , we will make the simplifying 
assumption that, experimentally, V A = V B = 1. Then, the total number of events is given by 

No - C i = C 

^£^p (167) 

assuming equal running in all four polarization modes. Hence, the number of events is cut in half, but the 
asymmetry is doubled. So, 

AAf, v 1 AAi , x 

Tjr-vor (168 » 

This means that we gain roughly a factor of V2 in precision by looking at Af[ V . One could argue, of course, 
that there are two Ai; A u i and Ai u . Then, one could combine them to improve the errors. Theoretically, the 
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two asymmetries are equal, but experimentally they only approach each other in the limit of infinite events. The 
problem is that these are not independent measurements. Hence, one cannot simply add the AAi in quadrature. 
If we define Ai as being the experimental average of the two Ai , a proper treatment of the errors yields 



AAj 1 
Ai x^Vo 

1 



l-A? V2 A PV - X 1 1 



+ 



APV2 ' APV 
A ll A ll 



-,1/2 

(169) 



In the limit of small A PV , AA x j A x -> AA PV /A PV but in general AAJAi > AA PV /A PV . Statistically, one 
therefore does best with A PV when two polarized beams are available. Although, Ai is not much worse for typical 
experiments, Z production in particular, where the fractional errors should be virtually identical in Ai and in 
A PV . Of course, the above analysis is only strictly valid in the limit of full beam polarization. Inclusion of 
partial beam polarization will not change our conclusions, however, since we can always 'pretend' that the beams 
are fully polarized. Then we are measuring the polarized parton distributions multiplied by the corresponding 
beam polarizations, rather than just the distributions themselves. When considering Xp (or rapidity) differential 
cross sections, Ai is the quantity of interest H] since it is more directly related to the ratio of the polarized to 
unpolarized parton distributions. 

Fig. |(b) gives A PV for 5 < M < 125 GeV. It has exactly the expected behaviour. Since AA PV ~ (1.5 ~ 
2%)/V, we can clearly distinguish between the possible solutions, and hence are quite sensitive to both Au and 
Ad in the region x ~ .18 . 

XI. W /± -BOSON PRODUCTION AT RHIC 

It is important to determine the intermediate- and large- a: behaviour of the antiquark distributions at high /i 2 
since this (along with the large- a; gluonic behaviour) influences the behaviour of the antiquark distributions in the 
limit x — > 1 at lower /x 2 , relevant to deep- inelastic scattering. This is because when one evolves from a low energy 
scale to a high energy scale, the large- a; behaviour influences the evolution at lower x, and vice-versa. Using W ± 
production at RHIC, we can gain insight into the x — > 1 behaviour of the antiquark distributions, where other 
experiments (such as deep-inelastic scattering) will have little or no sensitivity. This statement applies equally to 
the polarized and unpolarized antiquark distributions in the proton. 

W ± production at RHIC is ideally suited for this purpose because of the high event rate, resulting from RHIC's 
high luminosity, and because of the flavor specificity of the cross sections. In Figures ||(a) and ||(b), respectively, 
we plot the W + — > l + vi and W~ — » l~vi production cross sections for 200 < yfS < 700 GeV and for decay into 
one type of lepton (i.e. muons). Since one cannot measure M on an event by event basis, we have integrated over 
it. The HOC have the same structure as for 7* and Z production. At this point though, it is instructive to look 
at the variation of the QCD corrections with yfS. We understand the origin of the negative qg contribution, from 
the discussion in section IX. The fact that the relative magnitude of the qg contribution increases with yfS is 
understood to be a consequence of the increasing phase space in the inital state, reflected in a larger integration 
region over x a and Xb (or w). What we are really observing is the effect of varying r, since M ~ M\y for all y/S. 

One normally thinks of the large QCD corrections as arising mostly from the term ~ ir 2 in the qq subprocess, 
which comes from the virtual corrections. The rationale is that all the other terms are of order one and will tend 
to largely cancel among themselves. Clearly, things are somewhat more subtle than that though, since the term 



7 in (71) makes the finite virtual corrections small and negative. Also, the corrections to the qq subprocess 



clearly increase with decreasing y/S (i.e. increasing r). In order for the net correction to be positive, there mus t 



be a large, positive contribution arising from the qq bremsstrahlung. Such a large term is indeed present in ( 131 ) 
namely, the term ~ —4(1 + tu)m(l — w) which gives a large, positive contribution from the integration region 
near w = 1. There is also a smaller term ~ 41n 2 (l — wx)5(l — w) arising from the '+ '-distribution term. As 
t increases, these terms make a larger contribution relative to the other, potentially negative, bremsstrahlung 
contributions. The other term arising from the '+ '-distribution could give rise to sizable positive corrections as 
well, but its exact relative magnitude depends on the details of the parton distributions. In the above picture, 
the behaviour of the corrections is well understood. 
Considering the decay channel 

w ± -> M ± y , (no) 

for V# = 500 GeV, one predicts roughly 105,000 W + (fi + ) and 27,000 W~(n~) events. Not taking into account 
the details of the cuts, this corresponds roughly to an error on (double- and single-spin) asymmetries of 

AA W+ „ - 3% AA w- „, - Q % n?n 
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With such large rates at y/S = 500 GeV, it is not unreasonable to consider going to lower y/S as a way of probing 
larger x. In fact, the whole energy region, 

250 < VS < 700 GeV, (172) 

is interesting. Experimentally, the higher energies may be difficult to access. For the case of greatest experimental 
relevance (500 GeV), we are most sensitive to x ~ .16. 
The double-spin asymmetries adopt a simple form 

w + AuAd w - AuAd 

A ll — A ll — =1 ! \ LI ' S ) 

ud ud 

so that A^ + is sensitive to Ad and AW is sensitive to Am. From the signs of the polarized antiquark distributions 
for the respective sets, we expect 

GSA : < + < 0, AW~ > 0, 

GRSV : AW + > 0, AW~ < 0. (174) 

AW + and AW are shown in Figures ^|(a) and||(b), respectively. The signs are as expected. We also see that 
| -A/2 1 increases as decreases (i.e. x increases), which is a consequence of (150). This is a general trend which 



makes low energy measurements feasible. With high precision An measurements possible at yfS — 500 GeV, 
there is no problem in disentangling Am and Ad near x = .16. 

For AjY + , the HOC's have the same structure as A7, ' z , for the same reasons. For AW , the qg corrections 
have opposite sign since it is now the dg subprocess which enters. Another interesting feature of this and other 
asymmetries in production is that the effect of the HOC is comparable to the expected uncertainty in their 
measurement. Hence, a complete analysis must make use of the QCD corrections in order not to waste the good 
statistics. 

The W + single-spin asymmetry is roughly 

+ «AJ_ dAu = Ad _ An 
1 2nd 2ud 2d 2m v ' 

The second term is dominant over the first, but the first term does allow distinction between sets. We expect, 
based on the signs of the various Ad, 

|A^ + (GRSV)| > \AW + (GSA)\. (176) 

The W~ single spin asymmetry is approximately 

j^w~ dAu uAd Am Ad (177) 
1 ~~ 2du ~ 2du ~ ~2u~ ~ ~2d > ' 

Now the two terms are more comparable in magnitude, since u < d at large x and all the sets assume Am = Ad. 
Thus, the magnitudes of Am and Ad are important as is the sign of Am. In this case, we expect 

AW~ (GSA) > AW~ (GRSV), (178) 

since there is a cancellation in the GRSV case. In the GSA case, there is an enhancement rather than a cancellation 
and we expect large asymmetries at small y/S, where x is large. 

Figures 0(a) and^(b) present AW and AW , respectively. We see that the behaviour is exactly as expected. In 
all cases, distinction between the various sets is straightforward; more so for A^ , however. For AW + , knowledge 
of the HOC's are particularly important for this separation, especially at lower \/~S. For AW , we probe Am 
in a very direct fashion. The asymmetries are large throughout. They are roughly constant for the GRSV set, 
but increase in magnitude in the GS case as we go to lower energies, allowing a precise determination of Am 
over a wide range of x. Of course, no single measurement should be expected to determine any specific parton 
distribution exactly One must fit all the data, including the DIS data. 

Figures ||(a) and |](b) give A^ v for W + and W~ production, respectively. We notice again roughly a factor of 
2 enhancement over Ai, except near A^ = 1, where we are of course constrained by A^ < 1. Already, precision 
A\ measurements were possible, so A^ v , with the extra y/2 in precision, can pin down rather tightly the allowed 
sets of parton distributions. 

Z and W ± production at RHIC have been previously examined in LO as a tool for pinning down the polarized 
parton distributions in (ij^^]. In those studies (Z and VF ± ) rapidity differential asymmetries were considered, 
which arc quite useful in pinning down the x-dependence of the parton distributions. The general conclusions are 
the same. Also, transverse momentum distributions for Z and W ± production at RHIC were studied using Monte 
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Carlo methods which sum up certain bremsstrahlung graphs in |46[ . This observable is known to be sensitive to 
the polarized gluon distribution. 

Another interesting issue is the ratio u/d at large x and u 2 , which is unrelated to spin physics. Knowledge of 
the u/d ratio at large x and u 2 gives information on that ratio at even larger x at lower /i 2 , where there is very 
little experimental information. This was discussed at the beginning of this section. DIS is insensitive to the sea 
quarks at large x, since they are masked by the valence quark distributions. 

From the experimentally measured violation of the Gottfried sum rule , we can conclude that u(x) ^ d(x) 
for all x. Maximal violation of SU(2) flavor symmetry is usually taken to occur at larger x, where RHIC is sensitive. 
A typical assumption for x (and u 2 ) accessible to RHIC in W r± production is u/d ~ .5 jy]. This is based on 
the Pauli exclusion principle and explaining the Gottfried sum rule violation. The idea is that since the relative 
number of valence to sea quarks is increasing with x, the Pauli suppression effect will increase with x such that the 
u/d ratio decreases as there are more up valence quarks than down valence quarks, leading to greater suppression 
of the up sea. More recently, low energy fixed target Drell-Yan experiments at Fermilab have helped to disentangle 
d and u at lower /i 2 . 

The quantity of experimental interest (again considering the muonic decay channel) is 

_ a{W-) ^ du 

Rw = ^w+j * 3' (179) 

in crude approximation. Since the d and u contributions are quite well known, we directly probe u/d. With ~ 
130,000 events, we can measure Rw to high accuracy. At this level, we are again sensitive to the effect of the 
HOC's. 

Fig. |^ shows Rw for 200 < y^S < 700 GeV. Two sets are considered. First, the MRSG set and second, the 
MRSG set with d set equal to u. We see a clear separation between the two possibilities for all ^J~S. This will 
therefore be a crucial experiment for understanding the u/d ratio. This probe, at RHIC, has previously been 
studied in leading order and the same conclusions were drawn (see, for instance, [^4p9[). 



XII. THE FORWARD-BACKWARD DRELL-YAN ASYMMETRY 

As discussed in the introduction, the principal purpose of performing a high precision measurement of the 
forward-backward Drell-Yan asymmetry at Fermilab is to precisely determine sm 2 6w- In order to accomplish 
this, it is necessary to take into account the QCD corrections. We have seen that in the case of the spin dependent 
asymmetries in pp collisions, the qg subprocess often destabilizes the asymmetry. As we will see, for the forward- 
backward asymmetry, the corrections arise predominantly from the qq subprocess. The explanation for this 
feature will be given in detail. 

In Fig. |H](a), we present da/dM for l + l~ production in pp collisions at \/~S = 1.8 TeV. As usual, we include 
only one type of lepton-pair (i.e. u + fi~ or e + e~). Also, as in the previous numerical calculations, we use the 
MRSG [|l| unpolarized parton distributions in calculating physical cross sections. The peak in the cross section 
is quite a bit larger than at RHIC due to the higher energy (i.e. smaller x) and the valence-valence contributions. 
Here, and throughout, we use u = M . The effect of varying u will be discussed later in this section. 

Let us define the forward-backward (lepton) asymmetry as 



da^ F - B /dM 

AFB ~ d^/dW (180) 



where a F±B were defined in (|l2|). We expect large Afb in pp collisions in mass regions where Z — 7 interference 
dominates since the asymmetry arises, at the subprocess level, from axial - vector interference so that the Z 
contribution is pure axial and the asymmetry is unsuppressed. At the Z-pole, on the other hand, the Z vector 
couplings enter and lead to a suppression, so the asymmetry will be somewhat smaller there. We can see this 
explicitly by noting that the leading order subprocess level asymmetry is proportional to 

,qq 9al9aq9vl9vq f'ISI'l 

(9 2 al +9 2 Mq+9 2 V q) [ ' 

at the Z-pole. Hence there is an overall small factor, 

9vi = -i(l-4sin 2 M, (182) 

which is sensitive to sin 2 6w since sin 2 6w— -23, so that a small fractional change in sin 2 8w leads to a large 
fractional change in Afb- 

Fig. [l0](b) shows the effect of the QCD corrections on Afb- We see that, unlike the spin asymmetries (double- 
spin in particular), it is t he q q subprocess which accounts for the dominant corrections. We understand this in the 
following way. From Fig. [HJ(a), we see that the magnitude of the qg correction to a F+B is rather small compared 
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to that of the qq subprocess. Thus, in order for the qg subprocess to make an appreciable contribution to Afb, 
the corrections to o F ~ B would have to be quite different from those to <j F+B . Since one has the same parton 
distributions in both cases, all the difference arises from the differences in the subpro cess cross sections. As we 
have explained before, the dominant corrections come from the term ~ ln(l — w) in ( |132| ) which has the same 
form in both <j F+B and a F ~ B . Consequently, the qg corrections to a F amount to a multiplicative factor which 
basically cancels in the ratio, in the mass region of interest. We also expect small corrections to Afb arising from 
the qq subprocess, for the same reasons. The question then is: why are the qq corrections to Afb observable 
while the qg corrections are not? The correction basically arises from the differences in the hard bremsstrahlung 
contributions to a F±B , fc*|' F±B , given in (po|). We note that the difference vanishes in the limit w — > 1, relevant 
for the large M limit and for the dominant corrections to a F±B at lower M. This explains the vanishing of 
the corrections as M increases. At intermediate and low M, however, there will be some contribution to Afb 
arising from smaller w, in which case the term, k^ F ^ B = 2(1 + w)\nw, contributes negatively (oppositely to 
the Born term) in a significant fashion, thus accounting for the reduction in the magnitude of Afb- Now, for the 
qg subprocess, the differences in the k^ F±B all vanish in the limit w — ► 1 and are small at smaller w, thereby 
accounting for the relative smallness of the qg corrections to Afb- 

In the end, the corrections to Afb almost exactly amount to a multiplicative factor (less than I) which is given 
in Fig. [H]. We see that it increases with increasing mass, approaching unity. In the Z-po\e region, it is ~ .975 . 
This result is in good agreement with the finding of jl(J, considering that they use a different definition of Afb 
beyond LO, and a more involved approach, taking into account experimental cuts. Also, in their approach, one 
could not understand properly the structure of the QCD corrections, since the hard bremsstrahlung is handled 
via Monte Carlo integration. The spikes at M ~ 89.3 — 89.5 GeV simply reflect the fact that Apg° and Ap% 
vanish at slightly different points and intersect just above the zero. 

Figures |l2](a) and |l2|(b) show Afb in the Z-pole region for sin 2 9w= -2315. The effect of the QCD corrections is 
shown as is the effect of changing sin 2 9w by an amount of ±.0005. The magnitudes of the effects are comparable, 
but the character is distinctly different. Changing sin 2 9\v amounts to a shift in Afb rather than a multiplicative 
factor. As a result, changing sin 2 9w shifts somewhat the zero of Afb- Since the QCD corrections do not shift 
the zero appreciably, one sees that measuring precisely the zero in Afb may allow one to get a good handle 
on sin 2 9w without worrying about neglected higher order QCD effects (and possibly other uncertainties such 
as choice of parton distributions). How precisely this may be measured is left for a separate study. Using the 
GRV NLO (MS) set |5(J produced no appreciable change in the predictions. The CTEQ3M set |5^|, on the other 
hand, gave slightly different predictions for the forward-backward asymmetry. Hence, closer agreement between 
the various sets is required before precision determinations of sin 2 9w are possible. 

The fact that the zero of Afb does not change appreciably under HOC's is understandable since the zero 
essentially arises from the zero in the subprocess cross section, which in turn depends on M in a way independent 
of parton distribution effects. Each subprocess, at NLO, is proportional to a Born term with a well-defined zero. 
So any shift in the zero must arise predominantly from the fact that each quark flavor passes through zero at 
slightly different values of M . The NLO corrections however, do not appreciably affect the weighting of the 
various quark flavors relative to the LO weightings, since the same parton distributions enter. Hence there is very 
little shift in the position of the zero. 

We studied the effect of varying fi on Afb and it was found to be negligible. This is not unexpected since 
the /x-dependent part of the corrections has the same form in a F ~ B and a F+B . Hence, the major uncertainties 
are in the parton distributions, the neglected nonleading corrections, the QED corrections and possible intrinsic 
transverse momentum effects. Normally, one would not consider the latter effects at such high /i 2 , but as we are 
dealing with a very high precision measurement, they should not be taken for granted as being negligible. There is 
also a pure QED contribution to Afb at order a 3 , which was originally studied in connection with e + e _ — * 

With C = 110 pb , Afb for e + e pairs in the Z-pole region was measured with a statistical error of roughly 
±20% [ p3| . The maximum possible £ for Run II after several years running is 100 fb" 1 [0. We will take £ max = 
70 fb^ 1 as being a more realistic (if not optimistic) value in determining the best possible measurement of Afb- 
Then, statistically, we expect to be able to measure Afb to (at best) ±.8%. Hence, our statistical error goes 
down by a factor of roughly 25. The same statement applies to sin 2 9w- The statistical error on sin 2 9w of ±.003 
which is obtained from 110 pb -1 is reduced to ±.00012 with 70 fb _1 . If one takes into account two detectors and 
both muons and electrons, in a best case scenario, we could multiply our number of events by a factor of 4 and 
get ±.00006 as an error on sin 2 9w- For a more realistic 30 fb _1 , we get an error of ±.00009 ~ ±.0001. If we do 
not combine the data from both detectors, muons and electrons, this goes up to ±.0002. Either way, the error is 
very small. 

Of course, this error analysis is very naiive and all the systematic errors must be put under tight control for a 
realistic ultrahigh precision measurement as discussed above. A more detailed analysis is beyond the scope of this 
paper. Nonetheless, if any of the above scenarios could be realized, this would be the best sin 2 9w determination 
available. In fact, it would be better than the world average. Even with lower luminosity running such as 10 fb _1 
and not combining muons and electrons, the measurement is of the same precision as the SLAC measurement 
and hence could help to resolve the SLAC-LEP discrepancy. 
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XIII. CONCLUSIONS 



Complete analytical results for mass differential Drell-Yan type cross-sections relevant to all possible initial lon- 
gitudinal polarization states were calculated at the one-loop level in QCD. Interference between bosons of arbitrary 
mass, width and couplings was considered. For all observables considered, the corresponding forward-backward 
cross sections were determined. The results were presented in a form valid for all consistent n-dimensional regu- 
larization schemes. A survey of constraints on allowable factorization and regularization schemes was given. The 
mechanism behind scale dependences was discussed in some detail and in a general fashion. 

NLO predictions for all longitudinal Drell-Yan type processes at RHIC (W^, Z and 7*) were made using 
polarized parton distributions which fit the recent DIS data. The HOC's increased the cross sections substantially 
and had a major impact on the asymmetries, while preserving the features of the LO asymmetries. The exact sign 
and magnitude of the HOC's depended on the details of the polarized parton distributions used, especially the 
sea and gluon distributions. Faced with either low rates or small asymmetries, 7* production did not appear very 
interesting at face value, for longitudinal polarization. On the other hand, if the agreement between the various 
parton distributions at small x is accidental, rather than constrained, 7* production at RHIC could demonstrate 
that possibility by yielding an unpredicted asymmetry. The Z- asymmetries were all quite sensitive to the sea 
quarks; the parity violating ones being the largest, with unexpected sensitivity due to a coincidental cancellation 
between u and d valence contributions. With large rates and asymmetries, production can directly measure 
the polarized sea and valence distributions as well as the unpolarized u/d ratio. Lower energy running could 
directly measure u and d at rather large x (both polarized and unpolarized). 

NLO QCD predictions for the forward-backward lepton asymmetry at Fermilab were made. The QCD cor- 
rections amounted almost exactly to a multiplicative factor on the asymmetry which was a function of the mass 
of the lepton pair produced. This function was found to be less than unity throughout and approached 1 with 
increasing invariant mass. In the Z-pole region, the factor was roughly .975 . It was thus observed that the zero 
in the asymmetry was quite stable under QCD corrections. The zero was rather sensitive to sin 2 9w, however. 
This suggested an alternate method for determining sin 2 9w (he. measuring precisely the zero). The details were 
left for a separate study, however. Based on expected luminosities for Fermilab's Run II and previous Fermilab 
determinations which used a standard algorithm for extracting sin 2 9w from Aps, high precision sin 2 9w mea- 
surements were found to be possible at a level much better than (statistically) , or at worst competitive with, the 
best measurements presently available and should be able to resolve the SLAC-LEP discrepancy. 
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FIG. 1. (a) The cross section, da/dM, versus M, for 
production in pp collisions at ■s/S = 100 GeV; (b) 
corresponding double-spin asymmetry, An, for various sets 
of polarized parton distributions. Details given in text. 




FIG. 3. (a) The cross section, da/dM, versus M, for 
l + l~ production in pp collisions at yS = 500 GeV; (b) 
corresponding double-spin asymmetry, An, for various sets 
of polarized parton distributions. Lines as in Fig. [j]. 
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FIG. 5. (a) The total cross section for W + (— > l + vi) pro- 
duction, versus y/S, in pp collisions; (b) corresponding cross 
section for W~(-^> l~ v{) production. Lines as in Fig. hi. 
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FIG. 12. (a) The A FB of Fig. 0(b), but in the mass 
region, 80 < M < 100 GeV. The effect of varying sin 2 9 W 
by ±.0005, on Apg° , is shown, (b) Same as (a), but in the 
mass region, 88 < M < 92 GeV. 
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TABLE I. The AkTij which define the various factoriza- 
tion schemes of interest; here k — u,l. 
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